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SPECTRAL AND HODGE THEORY OF ‘WITT’ INCOMPLETE CUSP 

EDGE SPACES 

JESSE GELL-REDMAN AND JAN SWOBODA 


Abstract. Incomplete cusp edges model the behavior of the Weil-Petersson metric on 
the compactified Riemann moduli space near the interior of a divisor. Assuming such 
a space is Witt, we construct a fundamental solution to the heat equation, and using 
a precise description of its asymptotic behavior at the singular set, we prove that the 
Hodge-Laplacian on differential forms is essentially self-adjoint, with discrete spectrum 
satisfying Weyl asymptotics. We go on to prove bounds on the growth of L^-harmonic 
forms at the singular set and to prove a Hodge theorem, namely that the space of L^- 
harmonic forms is naturally isomorphic to the middle-perversity intersection cohomology. 
Moreover, we develop an asymptotic expansion for the heat trace near t = 0. 


1. Introduction 

On a compact manifold M with boundary dM which is the total space of a fiber bundle 

Z^dM-^Y, ( 1 . 1 ) 

with Z, Y closed manifolds, an incomplete eusp edge metrie pice is, roughly speaking, a 
smooth Riemannian metric on the interior of M which near the boundary takes the form 

gii,e = + x‘^^gz+ Tr*gY+g, k > 1, ( 1 . 2 ) 

where gy is a metric on the base Y, gz is positive definite restricted to the fibers, x is 
the distance to the boundary (to first order), and 5 is a higher order term. Thus near the 
boundary {M,gice) is a bundle of geometric horns over a smooth Riemannian manifold Y. 
When fc = 3, such metrics model the singular behavior of the Weil-Petersson metric on the 
moduli space of Riemann surfaces, as we discuss below. 

In this paper, we study the Hodge-Laplacian 

A := = d6 + 6d (1.3) 

acting on differential forms. Our first result shows that under conditions which contain the 
main examples of interest, one need not impose ‘ideal boundary conditions’ at dM in order 
to obtain a self-adjoint operator. 


Theorem 1. Let (M, g^ce) be an incomplete cusp edge manifold that is ‘Witt’, meaning that 
either dim Z = f is odd or 


Assume furthermore that gice 


= {0}. (1.4) 

satisfies (|2.9[) below and that the parameter k in (1.2) satisfies 


fc > 3. 


(1.5) 
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Then the Hodge-Laplacian acting on differential forms is essentially self-adjoint and 

has discrete spectrum. 


Thus, by the spectral theorem [35], there exists an orthonormal basis of Lf{pP{M)) of 


eigenforms 




We also prove that the distribution of eigenvalues satishes 


“Weyl asymptotics,” concretely, for fixed degree p 

#{j I < A^} = c„ Vol(M, 5 ,,e)A” + o(A") as A ^ oo. ( 1 . 6 ) 

See (4.2 for the proofs of Theorem]^ and of the asymptotic formula in (1.6). 

Having established these fundamental properties of the Hodge-Laplacian on such spaces, 
we turn to the next natural topic: Hodge Theory. Here the object of study is “Hodge 
cohomology”, or the space of harmonic forms, 


nl,{M,g,,e) = {a e L^{nP{M),g,,e) | da = 0 = Aa} , (1.7) 

and one phrasing of the Hodge theory problem is to hnd a parametrization for , gi^fj 

in terms of a topological invariant. As described in [H], in analogous settings the relevant 
topological space for Hodge theoretic statements is not the manifold M, but the stratified 
space X obtained by collapsing the fibration at the boundary over the base. 


X := M/{p ^ q\p,q& dM and 7 r(p) = 7 r(q)}. 


( 1 . 8 ) 


In (4.3 we will prove the following. 


Theorem 2. For a cusp edge space {M,gice) whose link Z satisfies the Witt condition (1.4), 
there is a natural isomorphism 


nl.{M,g,,e) - IH^{X), 


(1.9) 


where IHm is the middle perversity intersection cohomology of X. Furthermore, differential 
forms 7 S , gice) admit asymptotic expansions at the boundary of M. 

Moreover, if Z ^ S-^, the sphere of dimension f, then X is homeomorphic to a differen¬ 
tiable manifold and the isomorphism ( |1.9[ ) becomes 

nl2{M,g,,,)c^H*^j,{X), ( 1 . 10 ) 

where the latter is the de Rham cohomology of X. 


We recall the relevant facts about intersection cohomology, originally defined by Goresky 
and MacPherson in [mn], in Section |4.3| below. The equivalence in ( |1.9[ ) will follow using 
the arguments from Hunsicker and Rochon’s recent work m on iterated fibered cusp edge 
metrics (which are complete, non-compact Riemannian manifolds). To elaborate on the 
asymptotic expansion for L^-harmonic forms 7 , we will show in Lemma 19 below that in 
fact 


nl2{M,g,,,) = {a G L\nP{M),g,,,) \ A^-'a = 0} , 

(that the former is included into the latter is obvious), and we show that elements in the 
kernel of A®*^' have expansions at dM analogous to Taylor expansions but with non¬ 
integer powers, a statement which can be be interpreted as a sort of elliptic regularity at 
the boundary of M. 

One application of these results, and to putative further work we describe below, is to 
the analysis on the Riemann moduli spaces Ai-y^e. of Riemann surfaces of genus 7 > 0 with 
£ > 0 marked points. These spaces carry a natural Lf metric, the Weil-Petersson metric 
gwPi which near the interior of a divisor is an incomplete cusp edge metric with fc = 3. 
In general divisors may intersect with normal crossings, but in at least two cases only one 
divisor is present. 
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Theorem 3. Let (also known as the moduli space of elliptic curves) and A4o,4 be 

the spaces of, respectively, once punctured Riemann surfaces of genus 1 and 4 times punc¬ 
tured Riemann surfaces of genus zero, modulo conformal diffeomorphism. Then the Hodge- 
Laplacian /Spwp qji differential forms is essentially self-adjoint on Lf with core domain 
Cf°orb Section 12) with discrete spectrum and Weyl asymptotics, and if M.\p and A4o,4 
denote the Deligne-Mumford compactifications (see e.g. [HiiSB];. Then the de Rham coho¬ 
mology spaces HdR{Aiip) are naturally isomorphic to 'Hf 2 {M^i,i,gwp), o^nd the same holds 
for Mo,4- 

We discuss the proof at the end of Section though this is really a direct application 
of our results together with the recent work on the structure of the Weil-Petersson metric 
near a divisor in [2S] and |28) . 

This article is partly motivated by Ji, Mazzeo, Muller, and Vasy’s work [20j on the spectral 
theory of the (scalar) Laplace-Beltrami operator on the Riemann moduli spaces Mg, for 
which it was shown by methods different from ours that it is essentially self-adjoint and 
its eigenvalues satisfy a Weyl asymptotic formula. Here they analyze incomplete cusp edge 


spaces with normal crossings, and find in particular that the value A: = 3 in (1.51 is critical; 


indeed for values k < 3 one does not expect self-adjointness. It would be interesting (though 
more complex) to hnd a parametrization of the space of closed extensions of incomplete cusp 
edge Laplacians with k < 3, which is expected to be infinite dimensional, e.g. by [3]. 

In contrast with [^, since our eventual goal is Hodge and index theory on moduli space, 
our main technical contribution is the construction and detailed description of the heat 
kernel H = exp(—tA®*'’”). Indeed, our approach to establishing Theorem (which justifies 
the use of the word ‘the’ in the previous sentence) and Theorem]^ is to develop in Theorem 
|12| below a precise understanding of the behavior of a fundamental solution to the heat 
equation, which we only conclude is the heat kernel after using it to prove Theorem we 
establish asymptotic expressions for it at the singular set, uniformly down to time t = 0, 


obtaining in particular in Corollary 18 an asymptotic formula for its trace (which has 
potential applications to index theory, since our method for analyzing A®*'® may be used 
for other natural elliptic differential operators on these spaces as well) and fine mapping 
properties of A®*'”’ which allow us to analyze its kernel, i.e. harmonic forms. This is all 
described in detail in Section [H 

Essential self-adjointness of a differential operator P is typically a statement about the 
decay of sections u for which Pu € Lf. (Here the derivative is taken in the distributional 
sense.) The set of such sections is denoted 

2?max := I?max(A®‘") = {u & \ Pu & L'^] . (I.II) 

This is the largest subset of Lf which is a closed subspace in the graph norm ||u||r = 
||m||l 2 -f \\Pu\\l2. On the other hand, the smallest such closed extension from the domain 
Cf°{M) is the closure, i.e. the minimal domain 


,(A®-“) = \ u€ \ 3uk€ C^{M) with lim | 
I k—foo 


\uk - mIIt = 0 


(I.I2) 


The essential self-adjointness statement in Theoreml^says that the smallest closed extension 
is equal to the largest, i.e. that 

2?max = (1.13) 

and therefore there is exactly one closed extension. On the other hand, I?max is dual to I?min 


with respect to and thus if (1.131 holds then P with core domain Cf°{M) has exactly 


one closed extension, which we denote hy V = Vmin = 2^max and (P, I?) is a self-adjoint. 
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unbounded operator on 


Equation (1.13) is a statement about decay in the sense that 

decays fast enough near dM that 


to prove it we will show that a differential form a € T>n 
it can be approximated in the graph norm by compactly supported smooth forms. This we 
do using the heat kernel. 

Recall that the heat kernel H is a section of the form bundle 11: End (A) — )■ M° x M° x 
[0,oo), where M° is the interior of M and End(A) is the vector bundle whose fiber over 
{p,q,t) is End(A*(M); A*(M)), smooth on the interior M° x M° x [0,oo)t, which solves 

{dt + = 0 and Ht —Id, strongly as 11 0. (1-14) 

For a compactly supported smooth differential form a, the differential form 


/3(a;,t):= / H{'w,'w,t)a{w) dVolg 

JM 


■\w) 


solves the heat equation {dt + = 0 with initial data /3|t=o = a. One consequence of 

our precise description of H in Theorem below will be the following. 

Theorem 4. On a Witt incomplete cusp edge space {M,gice) with metric satisfying the 


assumptions in (2.9) below together with (1.5), there exists a fundamental solution to the 
heat equation Ht = H(w, w,t) in the sense of { 1.14[ ) such that for t > 0 

Ht : L2(M; n*{M)) (1-15) 

and such that Ht and dtHt are bounded, self-adjoint operators on L^. 

Theorem 1^ implies the essential self-adjointness statement; indeed the fundamental solu¬ 
tion Ht directly gives a sequence (indeed a path) of sections on I?niin which approaches a 
given form in I?max- Namely, 

a G T^max Htct —>■ a in I?niin as 11 0. (1-16) 

As we see now, the proof of this is straightforward functional analysis given the conclusions 
of Theorem |4l 

Proof of essential self-adjointness using Theorem^ The proof has nothing to do with the 
fine structure of incomplete cusp edge spaces, it depends only on the soft properties of the 
fundamental solution H in Theorem To emphasize this, let (M, g) be any Riemannian 
manifold and P a differential operator of order 2 acting on sections of a vector bundle E 
with hermitian metric G, such that P is symmetric on L‘^{M;E). For t > 0, let Ht be a 
smooth section of End(E) M x M which depends smoothly on t and satisfies 


{dt P)Ht = 0 and 


lim Ht = Id, 

t-fO 


(1.17) 


where the latter limit holds in the strong topology on , and furthermore such that Ht 
and dtHt are self-adjoint on L^. 

Let u € VniaxiP), i-e. u G Pu G L^. We will show that u G T’min)^’) as well, and thus 
Indeed, we will show that 

-2 


HtU —>■ u in I?maxj he. that HtU —>■ u and PHtU —>■ Pu in H 


(1.18) 

This suffices to prove that u G I^min since HtU G I^min by assumption and I?min is a closed 
subspace of I?max in the graph norm. To prove (1.18), we note first that HtU —>■ it in 
trivially since Ht —>■ Id in the strong topology on Lf'. Also note that since u G 'Dmax, 
Pu G so HtPu —)■ Pu in also. Of course, this is not what we want; we want 
PHtU —>■ Pu, but in fact we claim that 


M G 


PHtU = HtPu, 


(1.19) 
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which will establish (1.13). 


It remains to prove (1.19). Note that for u £ (C>max and v £ L 
{Pu, Htv)L 2 by self-adjointness of Ht on while {Pu, Htv) l 2 = {u 


then {HtPu^v) 1,2 = 
PHtv)L 2 . Indeed, the 


adjoint domain of is 2?max, so for any / e G 2?max, {Pf,g)L^ = {f,P9)L^- But, 

then since PHt = —dtHt we see that {HtPu,v) 1^2 = —{u,dtHtv)i^ 2 . But dtHt is self-adjoint 
onL^ so {u,dtHtv)L 2 = {dtH tU,v) L 2 = -{PHtU,v)L 2 , and thus {HtPu,v)L 2 = {PHtU,v)L 2 
for all u £ (Dmaxjii € , i.e. ( 1.19|) holds. □ 


The central vehicle for the construction of the heat kernel is the construction of a manifold 
with corners via iterated radial blow up of the natural domain of the heat kernel, 

namely the space M x M x [0, oo)t; thus the interiors of these two spaces are diffeomorphic, 
and the blowup process furnishes a ‘blowdown’ map 

/?: ^ M X M X [0, oo)*, (1.20) 

which encodes deeper information about the relationship between the various boundary 
hypersurfaces (codimension one boundary faces) of and those of M x M x [0, oo)*. 

The upshot is that the heat kernel if, which lives a priori on the latter space, pulls back 
via (3 to be “nice” (precisely to be poly homogeneous, see Appendix [A|) on ^heaf In fact, 
in Section]^ we will construct a parametrix K for the heat equation directly on To 

obtain the actual heat kernel H we use a Neumann series argument to iterate away the 
error. 

The latter process builds on what is now a substantial body of work on analysis (in 
particular the structure of heat kernels) on singular and non-compact Riemannian spaces, 
going back at least to the work of Cheeger on manifolds with conical singularities [ziisiin]. 
Our approach here is more closely related to Melrose’s geometric microlocal analysis on 
asymptotically cylindrical manifolds m (a non-compact example) and Mooers’ paper |32] 
on manifolds with conical singularities (an incomplete, singular example). The general 
procedure, which one sees in both the parabolic and elliptic settings, is to express the 
relevant differential operator as an element in the universal enveloping algebra of a Lie 
algebra of vector fields, and to ‘resolve’ this Lie algebra via radial blowup of the underlying 
space. 

Two papers closely related to our work are, first, Mazzeo-Vertman m, in which the 
authors study analytic torsion on incomplete edge spaces, which are the k = 1 case of 
incomplete cusp edges. (They analyze the behavior of elliptic operators which have the 
same type of degeneracy in their principal symbols near the singular set as the Laplacian, 
so a welder class than considered here.) Their analysis involves a heat kernel construction 
using blowup analysis, which is slightly simpler in their context as the resolved double space 
has one less blow up (and thus the triple space is much simpler). One substantial difference 
between the incomplete edge and incomplete cusp edge cases is that in the incomplete edge 
case the space of self-adjoint extensions is substantially more subtle. In particular, a Witt 
space (this is a topological condition and has nothing to do with the value of k) that is 
incomplete edge may have infinitely many self-adjoint extensions if the family of induced 
operators on the fibers have small non-zero eigenvalues [3]; on incomplete cusp edge spaces 
the small non-zero fiber eigenvalues do not contribute. On the other hand, one expects that 
the zero mode in the fiber (the space of fiber harmonic forms) makes a simlar contribution 
in both the cusp and cone cases, in particular that an incomplete cusp edge space which 
is not Witt will have an infinite dimensional space of closed extensions on which ‘Cheeger 
ideal boundary conditions’ must be imposed to make the operator self-adjoint, as is the case 
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A second closely related work is Grieser-Hunsicker m, which uses also inhomogeneous 
radial blowups, in this case to construct a Green’s function for elliptic operators on a certain 
class of complete Riemannian manifolds (called ‘(/(-manifolds’) which require similar analysis. 
There are many other related works in a similar vein including, just to name a few, Albin- 
Rochon [3], Briining-Seeley [B], Gil-Krainer-Mendoza m, Lesch [m, Schultze [M], and 
Grieser’s notes on parametrix constructions for heat kernels M- For analysis of moduli 
space, to give just a sample recent work, we refer the reader to the papers of Liu-Sun-Yao, 
for example [SUES]. 
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2. Incomplete cusp edge differential geometry 

As described above we work on a smooth manifold with boundary M whose boundary is 
the total space of a fiber bundle with base Y and typical hber Z. (Again, we assume that 
Y, Z are both compact and without boundary.) We fix once and for all a boundary defining 
function x, meaning a function x € with {a; = 0} = dM and dx non-vanishing on 

DM, and thereby fix (after possibly scaling a; by a constant) a tubular neighborhood of the 
boundary, 

Uc:idM X [ 0 , 1 )^. (2.1) 

Note that near each point q on the boundary, there is a neighborhood V in dM oiy = TT{q) (tt 
is the projection onto the base) such that 7 r“^(F) ~ Fx Z, where ~ denotes diffeomorphism. 
Below, we will say that we work ‘locally over the base’ when we fix our attention on a 
neighborhood of <7 in M of the form 7 r“^(F) x [0, l)^;. If one chooses local coordinates y on 
the base and z on 7 r“^(g) ~ Z, then 

{x, y, z) form a coordinate chart on M in a neighborhood of q. 


( 2 . 2 ) 
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Let 


f:=dimZ, 6:=dimy. 


(2.3) 


In the space X defined by collapsing the fibers Z over the boundary (see (1.81), the image 
of dM via the projection tt: M —>■ X defines a subset 7r(9M), which can be identified with 
y, and TT{dM) forms the unique singular stratum of the stratified space X. (See, for instance 
[2] or [in], among many others, for a discussion of stratified spaces.) In the case Z ~ S-^ a 
sphere, X is homeomorphic to a manifold. In general, given p € 7r(dM), there is a basis of 
open sets at p consisting of canonical neighborhoods homeomorphic to 

V X Cs(Z), (2.4) 

where V is an open subset of V containing p and 

Cs(z) = [0,(5) X y/{o} X y. 

2.1. Differential form bundles and the Hodge-Laplacian. We will consider differential 
forms and vector fields which are of approximately unit size with respect to Riemannian 


metrics of the type in (1.2). These are the incomplete cusp edge forms, which are sections 
of the incomplete cusp edge form bundle, ‘'^®A*(M), whose smooth sections are generated 
locally over the base by the forms 

dx, dpi (i = 1,..., 6 = dimy), x'^dza (a = 1, ...,/ = dim 
Locally over the base, we have the isomorphism 

‘“A*(M) = A*(y) A x^^h*{Z) ®dx^ X*{Y) A x'^^K*{Z), 

where N is a multiplication operator on h*[Z) which takes a form a of degree q to qa. 
Note that with respect to this decomposition, if N is the operator acting on A'^(y) by 
multipliction by N, the exterior derivative d acts as 

0 


(2.5) 


( 2 . 6 ) 


d = 


d^ + {-lY^x-^d^ 




(2.7) 


Moreover, we will use the space of vector fields which are locally C°°{M) linear combi¬ 
nations of the vector fields 


dx, dy 


„—k 


d. 


( 2 . 8 ) 

=r*(M) = 


These vector fields are local sections of a bundle *“T(M) which is dual to 
We denote sections of ®“T(M) by Vice- 

We consider metrics g on M which are sections of Sym°’^(*'^®T*(M)), and which locally 
over the boundary take the form 

10 0 \ \ / dx 

g={dx dy’' x’’dz°‘ ) \ 0 (hij) 0 | + 0(a;'",go) | ( dy’ j, (2.9) 



0 (ka^) 


x^dz° 


with hij and and independent of x. Here O(a;^,go) refers to a 0{x^) norm bound with 


respect to an exact incomplete cusp edge metric go as in (2.10) below, and furthermore 
we assume that the O{x^,go) term is polyhomogeneous conormal, a regularity assumption 
defined precisely in Sectionj^below, which roughly speaking means that the coefficients have 
an asymptotic expansion at cc = 0 analogous to a Taylor expansion but with non-integer 
powers and with precise derivative bounds on the error terms. Metrics satisfying these 
assumptions are what we refer to henceforth as incomplete cusp edge metrics. (Note 
that the assumptions on g are stronger than merely assuming that g G Sym°’^(*“T*(M)), 
as the latter space contains e.g. x{x^dz ®sym dx), which does not obey the error bound.) 
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Remark 5. As is shown in |26j (see the introduction for further discussion) with previous 
results for example in [Ml ED the Weil-Petersson metric on moduli space takes the form 
(2.9) near the interior of a divisor and satisfies the polyhomogeneity assumption. 


As an example of such a metric, one can take an exact incomplete cusp edge metric, 
constructed as follows. Let hi S Sym°’^(y) be a metric on the base, and h = tt* hi its x- 
independent pullback via the projection tt : dM —> Y extended to our tubular neighborhood 
U ~ dM X [0, l)a; in the obvious way. The tensor k we take to be the pullback to of a 
tensor on dM which by abuse of notation we write as k G SymP’^{dM), having the property 
that k restricted to any fiber is positive definite. The metric = dx^ + h + cc^^k is an exact 
incomplete cusp edge metric. Note that in coordinates (x, y, z) such a metric takes the form 


go = { dx dy'^ x^dz^ 



( 2 . 10 ) 


Thus an incomplete cusp edge metric diffes from an exact one by a polyhomogeneous error 
of order 0{x^) in norm and can therefore be taken to correspond uniquely to an exact cusp 
edge metric. 

Note that an exact incomplete cusp edge metric (and thus an incomplete cusp metric) 
gives rise, locally over the boundary to well defined notions of the operators dy, dz, the 
Hodge star operators and *z, obtained by allowing the Hodge star associated 

with {Y,h), the Riemannian structure on the base F, to act on the A*{Y) factors in the 
decomposition (2.6) and *(z.k); the j/-dependent Riemannian structure induced on the fibers, 
on the x^^A*{Z) factors. (So, for example, *zx^dza = x^d-Ak dza.) Thus, we may also 
define unambiguously (still locally over the boundary) the Hodge-de Rham operators and 
dz- Indeed, a homogeneous differential form a, in the coordinates {x^y,z) can be written 
a = {f{x,y,z)dy^) A x^Pdz^{Adx) = ±{f{x,y,z)x'^Pdz^) A dy^{Adx), where / and A are 
multi-indices and |/|, resp. |A|, is their order. One can act via By = dy + (—l)^(l^l“i)+i 
dy-ky on f{x,y,z)dy^, where -ky is the Hodge star operator of {Y,h), and thus define 
or act via Bz = dz + (—l)t'(l^l-i)+i *2 on f{x, y, z)x^Pdz^ where *2 is the Hodge star 
operator of {Z, k), and thus define dza. Note that dz is thus a family of operators on the 
fiber Z parametrized by the base point y. The Hodge-de Rham operator can be decomposed 
locally over the boundary in terms of the base and fiber Hodge-de Rham operators according 
to the following proposition, which we prove in Section |2.2| below. 


Proposition 6. Locally over the boundary (see below (2.1)^, the Hodge-de Rham operator 
decomposes as 


= do + 5zPi + P25z + 


( 2 . 11 ) 


where do acts on sections of decomposed as in (|2.6|) by 


JO = 


-knZ I piY _Q _ 


-f i 
kx-^N 


(/-N) 
- d^ 


d^ depends on the base Y parametrically, acts in only the fiber direction and only on the 
A*{Z) factor, and is equal to the Hodge-de Rham operator for the Riemannian manifold k\y, 
and where the Pi are polyhomogeneous, bounded endomorphisms on *“A*, and E is as in 


(2.12)-(2.14). 
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To describe the term E, we say first it is a differential operator of order one which does 
not increase the order of blow up of polyhomogeneous distributions; in particular 

£;(a:S) = 0{x^), (2.12) 

for 7 a smooth, bounded section of “®A*. In fact, it is a b-differential operators on ice-forms 
with polyhomogeneous coefficients 

(2.13) 

Letting Vb,phg denote the polyhomogeneous vector helds tangent to the boundary dM, (2.13) 
means that E lies in the algebra of differential operators generated Vb,phg- Concretely, it 
satisfies 

E = axdx + Edyi + c°‘dz<^ + d, (2-14) 

for polyhomogeneous, bounded endomorphisms a,E, 0 °^, d, and where repeated indices are 
summed over. In general, an element Q G Diff((*pj^g(M;*®®A*) also satisfies (2.12), and is 
given locally by polyhomogeneous linear combinations of xdx,dy,dz, i.e. 

Q = X! ai^a,0{xdxydyd^, 

i+\a\ + \l3\<m 

where ai^a,p is a polyhomogeneous bounded endomorphism of *®®A*. 

The Hodge-Laplacian A = = d6 + Sd can now be decomposed along the same lines. 


Proposition 7. Locally over the base, A can be decomposed as follows 

A = Aq -f x~^P + x~^E, 


where Aq acts on forms decomposed as in (|2.6|) by 

Aq = Id2x2 


-dl - ’^dx + 4 ^ a ^’^ 




/fcN(l-fc(/-N))x-2 

-2kx-'^-^S^ 


A' 

-2kx-^-^d^ 
k{f -N){l-kN)x-^ 


(2.15) 


(2.16) 


where 


P = 6zQi + dzQ 2 + Qd.5z + Qidz 

and where the Qi G Diffj', pjjg(M) and E G Diff^ pj^g(M), as described below Proposifion 

2.2. Proof of Propositions!^ andTo prove Proposition!^ we will compare the Hodge- 
de Rham operator for g to that of the locally defined warped product metric 


Pp = { dx dy^ x^dz° 



(2.17) 


We compute the Hodge-de Rham operator for pp now. First of all, the Hodge star operator 
of pp, -kp, acts on forms decomposed as in (2.6) as 


(2.18) 


0 (-l)’^(^-9^)*y*2 

where qy denotes multiplication by the degree in Y and N denotes multiplication by the 
degree in Z, and b = dimF. Using ( !2.7| ), we see that the Hodge-de Rham operator for pp 
acting on fc-forms. Bp = d+ (—l)^l”“C*p d*p satisfies (2.11). We now use this to prove the 
proposition. 
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Proof of Proposition^ A brief calculation shows that (2.11) holds with Pi = P 2 = E = 0 
for the warped product metric Pp, defined locally on M in (2.17). Next, we claim that the 
general case follows from that of a warped product metric Pp. To see this, it suffices to show 
that if -kp continues to denote the Hodge star operator for pp, then the pointwise norm 


< Cx\ 


(2.19) 


Indeed, if this holds, then writing 


± (5 = * {dx A dx (B dy © dz) *, (2.20) 

we consider 6 — 6g^, the latter being the dual of d with respect to dg^, which is given by ±1 
times 


(kdx A dx k — kp dx A dxkp) © (kdy k — kp dykp'j 0 (kdz k — kp dzkp). 


( 2 . 21 ) 


The operator dz, as an operator on sections of satisfies dz = where 0 are 

endomorphisms satisfying m = 0{x~^). (Proof: dza = '^^.{x^dza) A x~^dz^a and x^dzi 
is a bounded endomorphism since x^dzi has unit length.) By the same rationale, dx A dx 
and dy are differential operators with bounded coefficients. By (2.19) it follows that in the 
difference 


5 — 6p = i(*c? * — d*p) 

= zL6d{* — *p) i (* — *p)d(* — *p) i (* — *p)d * 
= ±(5 * (* — *p) ± (* — *p)d{* — *p) ± (* — *p) * i5 


the terms coming from dy and dx derivatives have coefficients bounded to order x^, while 
those coming from dz have coefficients with bounded derivatives. Hence (2.11) follows. 

It remains to prove (2.19). First note that 

g®k-(»l — g^. . . (g) . 

and similarly for pp, then 


under the assumptions above, if we let 
p~^ denote the metric induced by p on 0 


9 




9: 


P 


< Cx'^ 


( 2 . 22 ) 


where the norm is the one induced by p on the tensor powers. Indeed this is obvious for 

^), and (2.22) follows. Now let 

I = {ii, 


I = 0. Moreover, p = pp + 0{x^) implies p ^ = Pp^ + 0{x 
ti be an oriented orthonormal basis of 1-forms for p, and for 


,ik) denoting an 

ordered subset of {1,... ,n}, let e/ = A • • • A be the corresponding fc-form. Then 
kei = e/c where E denotes the complement of I in {1,..., n} such that e/ A e/c is positive. 
Thus, we want to show that 11*6/ — kpCiW^ < For two ordered multi-indices I, J, letting 
Sj j = 1 ii I = J and zero otherwise, consider 


{kpCp,e-j)gp — { {kpCj,e.j)g^ dVolpp, dVolg^) g^ 

= (e/Aej,dVolgJg^ 

= ±(5/c,j(dV0lg,dV0lgJgp 

= ±<5/=.j(l + 0{x^)), 


(2.23) 


so since (*„ 


e j}g = I 

bound in (2.19|) holds. 


e/, ej)gp + 0{x^) = ±Sioj{l + 0{x'^)) 


{kei,ej)g + 0{x'^), the 

□ 


To see that Proposition!^ holds, simply use A = together with Proposition]^ 
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2.3. Fiber harmonic forms. We now discuss forms which are approximately harmonic 


with respect to the family of metrics k induced on the fibers Z by an ice-metric as in (2.91. 
We begin by working directly on the fiber bundle tt : dM — > Y with a submersion metric, 
i.e. a metric = 7r*/ii -|- k (so hi is a Riemannian metric on Y and k is assumed only to 
be non-degenerate when restricted to the fibers). Then a G fi'^(dM) is fiber harmonic 
if is harmonic with respect to k| 7 r-i(p) for every p € Y, i.e. if it lies in the kernel 

of One can obtain fiber harmonic forms as follow. Let a G r(Y; where 

is the bundle whose fiber over a point p G Y is the space of q-forms on 7r~^(p). Since a 
metric is present, there is a natural way to extend such a 5 to a differential form on DM; 
consider the space := L'^['k* 0.‘^{Y)) C L'^{dM; A‘‘) which is the closure of the space 
of differential forms pulled back from Y, and note that there is a unique differential form 
on dM which is at every point orthogonal to H'^ and which pulls back to a on each fiber. 
(This is analogous to being a horizontal vector field; vertical vector fields are defined in 
differential terms and the horizontal vector fields are defined with the metric. Indeed, being 
a vertical form here means being dual via the metric to a vertical vector field). Moreover, 
if a is harmonic on each fiber, we claim that a := ay is fiber harmonic. Namely, choosing 
a local trivialization of dM near p gY and coordinates {y,z), for a family of q-forms a on 
the fibers, the vertical extension satisfies 


ay = 


\A\=q 


dzA 


E ^ 

ld + |S|=9.|/|>l 


I,B 


dyi A dzB, 


(2.24) 


with a = X]|A|=g cL^dzA- Choosing another local trivialization(?/, z) we have z = z{y,z), y = 
y, dzi = {dz^/dz^)dP + /’'dyi. Thus the pullback of ay is the change of variables for 
differential forms on the fiber, on which A'^’^ acts invariantly, proving the claim. 

Let ker(A'^) —>■ Y denote the bundle whose fiber above each y is the kernel of A'^“’’'«. 
This bundle has a grading by the form degree, 

nl := ker(A^ : fl^Z) —^ fl«(^)), 

and by the above paragraph the vertical extensions of these forms also forms a vector 
bundle 'H'gy C A'^{dM), the sections of which are fiber harmonic. We now define the 
approximately fiber harmonic differential forms H to be the direct sum of the spaces 

n = where W := tt* A{Y) A (2.25) 

or, in words, the sections of form the space of differential forms on ~ dM x [0,1)^; 
pulled back from Y wedged with the vertical extentions of harmonic q-forms on the fibers, 
thought of as living on U and weighted by (so as to make them unit size.) 

The main point to note about sections of % in ( |2.25 1 is that a form ^ gH satisfies that, 
locally over the base, we have 

= 0{x^), (2.26) 

where A^’^ is the operator defined (again only locally over the base) in Proposition [?{ and 
as usual the 0{x^) bound is a pointwise norm bound. Indeed, this follows since by (2.241, 
under a change of local trivialization of the boundary fibration, we have 

x^dz, = x^{dzydV)dzy Y 0[x'^), (2.27) 

so x^y times 7 is 0{x^) times the g-form part of the pullback to the fiber, together with the 
fact that the g-form part of the pullback is harmonic. 

Below we will often work with forms which are merely polyhomogeneous (as opposed to 
smooth). These are sections 7 of Mphg(‘“A*), defined rigorously in (A.3) below. As described 
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above, these have an expansion at a; = 0 analogous to a Taylor expansion, but with non- 
interger and possibly also log^(a:) factors. In particular we will be forced by possibly only 
polyhomogeneous regularity of the metric to work in the larger space defined 

to be the subset of ^phg('“A*) such that 

x*Mphg(H«) C W) 0 (2.28) 

This space contains, in particular, sections of and also ice-forms 7 which are poly¬ 

homogeneous and can be written as 7 = a;'’°7o + 0{x^°~^^), where 70 S but moreover 
contains x^°^i + for 7 ^ € and indeed any polyhomogeneous ice- 

form whose expansion terms not bounded in norm by a:®o-i-fc have fiber harmonic coefficients. 
Thus, 7 G x^°Aphgi'H'^) implies that locally over the base, 

A^’«7 = 0(a:"'>+'=). (2.29) 


It is in fact possible to define an operator A.fib on certain forms which is an invariant 
version of the locally defined for example on smooth sections of ''^®A*(M), if one 

notices that a section a G defines, by restriction to the boundary, a 

smooth section of A*Y<S)^*fib — Y- Indeed, this follows since (2.27) and local computation 


produces a form on Y with values in fl*{Z). Letting a denote this section, one can then 
define a to be /3i = a:^'^(A^’^ a)v, the vertical extension of the fiber-wise Laplacian applied 
to the n*(Z) factor of 5, and weighted by to make it a smooth section of '™A*(M). 
Then /3i — = 0(x) locally over the base. One can in fact iterate this to find a section 

Afiba := /? such that /3 — A^’^a = 0{x^), since if xa' is a smooth section then letting a’ 
denote its value (as an 0*(Z)-valued form on Y) the ice-form j32 = x ■ x^^(A'^’*'5')v, then 
(/3i + /32) — A^’^of = O(x^), and this can be iterated up to the error O(x^) (and since all 
steps introduce O(x^) errors, no better.) We will let Afiba = j3. 

We will also need to ask when a, say, smooth ice-form can be written as A fib of another 
form, or equivalently, when, given a ice-form a, there exists a “^®A*-form 7 such that locally 
over the base a — A^’^j = O(x^). A related question is whether there is an invariant 
definition of projection onto the space of fiber harmonic forms. Indeed, given a a smooth 
section of ‘'^®A*(M), again considering the valued form on Y, a, over the boundary, if 
n*' denotes the projection onto the kernel of A^’'^, we can consider x*'^(n^a), and again 
can iterate this fc-steps down the Taylor expansion of a to obtain a projection 11^ onto the 
fiber harmonic forms. In particular, if a G x‘’°'Hpb,g, the definition of 11^ can be extended 
to 

Un : ^ + 0{x'‘), ( 2 . 30 ) 

where 11 ^'^ is the projection onto ker(A'^’^) defined locally over the base, so that, again 
locally over the base we have 


e x'=Diff^(M), 

and thus since A^’*^ = d^S^ + S^d^, 

fjHAZ,VYiH g 


(2.31) 


(2.32) 


where Diff^ is defined in (2.13). Moreover, it is straightforward to show that if a is a smooth 
section of x®°(‘“A*(M)) then 

n^a = 0(x"“+'=) ^ 37 such that A^’*' 7 -a = 0(x"“+'=), 


this holding locally over the base, with the form 7 defined globally. 


(2.33) 
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3. The heat kernel 


In this section we construct a manifold with corners as in (1.20) together with 

a fundamental solution to the heat equation which is a polyhomogeneous conormal distri¬ 
bution on with prescribed leading order terms in its asymptotic expansions at the 

various faces (see Theorem 12). To do so, after the construction of we perform a 

parametrix construction and then use this parametrix to obtain the fundamental solution 
itself via a Neumann series. 


3.1. Heat double space. The space is obtained by performing three consecutive 

inhomogeneous radial blowups of M x M x [0, oo)t. Such blowups create, from a given 
manifold with corners X and some other data including a submanifold N, another manifold 
with corners [X\ N]inhom, which is diffeomorphic to the complement of N+{N). Here N+{N) 
denotes the inward-pointing normal bundle of N which we think of as an open neighborhood 
of in X; thus [X] N]inhom is a manifold with corners with one more boundary hypersurface 
(bhs) than X. It comes together with a blowdown map /3: \X]N]inhom. —^ X which is 
a diffeomorphism of the interiors and satisfies that C°°{[X]N]inhom) contains (3*C°°{X) 
properly, i.e. there are functions which are not smooth on X but nonetheless pull back to 
smooth functions on the blown up space. 

In more detail, the submanifold is a ‘p-submanifold,’ meaning that there are local 
coordinates {x, x', y, y') € x with x € y G in which N is defined 

(locally) by X = 0, y = 0. The space [X, N]inhom with homogeneities xq ~ x“^ ^ ... ^ Xp” ^ 
~ ~ y^’ with ttj, /3fc G N, 1 < ai < • • • < cip, Pi <■■■< j3q is a manifold with corners 

whose set of boundary hypersurfaces contains that of X naturally, and [X, Njinhom das one 
new boundary hypersurface, which we call nf for ‘new face’. Assuming for the moment that 
ctp > Pq and also that ap\aj,ap\Pk for all j,k, the space PA([X, N]inhom) is by definition 
the set X \ N U T where T is the set of paths 7 ( 5 ) in X with 7 ( 0 ) G N satisfying that 
Xp = ap{s)s for smooth non-vanishing s, and Xj = aj{s)s°‘’‘/°‘Pyk = bk{s)s°‘p^^'‘ for smooth 
aj,bk, all other coordinates being smooth functions of s, modulo the equivalence relation 
7 i ~ 72 if tde coordinate functions agree to order higher then the stated vanishing order 
(e.g. Xj( 7 i(s)) — Xj(j 2 (s)) = The space X\A^UF is naturally isomorphic to 

X \ A^U A^“''(A^) and can be given a smooth structure so that the following polar coordinates 
are smooth and valid on a collar neighborhood of the introduced face nf, 

Pnf = (xo+x^^ +---+x^’’ + \yif^ -b •• • -b \yqf'>)^^°‘’’, 


Pnf = ( 



xi yq , 

Op/ai ’ ■ ■ • ’ OLp/P, > 

Cnf Pnf 


A full set of (polar) coordinates is then (pnf, ?!>nf, yO 

For a detailed definition of such spaces we refer to Melrose’s work |29l Chapter 5] which 
contains a more general construction which does not assume that one has in particular a 
fixed extension for the manifold N away from the boundary, (whereas here we fix once and 
for all a boundary defining function x as in ( 2 . 1 ), which will give all the desired extensions 
below). See also [13HI]. 

(1) We first blow up the fiber diagonal in the corner. This is the subset of dM x dM x 
{0} C M X M X [0, oo)t consisting of points (p, q, 0) with 7 r(p) = 7 r(y) where tt is the 
projection of the fibration dM onto its base. If local coordinates (x, y, z) are chosen 
as above, this set can be written 


Bo := {x = x = t = disthiy, y) = 0} , 


( 3 . 1 ) 
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where dist;i(*, •') is the distance function on the base (Y, h). In fact, Bq is naturally 
isomorphic to dia.gf^y^{dM) x {0} where diag£|2(9M) is the “fiber diagonal,” i.e. the 
fiber product of dM Xgb 9M. We let 

^heat,l := [M X M X [0,00)t;Bo]inhom, (3.2) 

with t \y — y\^- To be precise, ^ is the parabolic blowup in time 

of the set Bq as defined in EH Chapter 7]. In particular there is a blowdown map 
Pi: ^ X [0, oo)t, and polar coordinates on ^ near Pp^{Bo) (once 

coordinates y, z are chosen on dM) are given by 

p= (^t + x"^ +P + \y- yp), 

(3.3) 

\p^ p p p J 

= {Pt,Px,Ps,Py), along with y, z, z. 


The set {p = 0} is a boundary hypersurface on ^ introduced by the blow 

up; we call it ffi; we will see that only the projection of the heat kernel onto the 
zero mode in Z is relevant at the face ffi. Letting s = x/x, the interior of ffi 
is the total space of a fiber bundle over Y x (0,oo)s, which is the fiber product 
dM Xfib dM Xflb TY X Kj/ where t' is a rescaled time variable (see ( |3.5| ) below). 
Indeed, the map from ffi to the base Y is simply Ppa^ 

(2) The preceding blow up does not resolve the term arising from (2.16). To 

accomplish this, we blow up the subset of ffi defined in polar coordinates by 


Bi := {p = = (j)y = 0, px = Ps} , (3.4) 

i.e. hy p = 0,(j) = (0, l/-\/2, l/-\/2, 0), inhomogeneously so that near the new face, 
ff, the function is smooth, and furthermore so that tdp. is non-degenerate, the 

latter condition being satisfied if {x — x)/y/i is smooth up to the interior of ff. Near 
Bi we can use projective coordinates 


X, s = xjx, p = ^ , t'=tjx‘^, (3.5) 

X 

along with y, z, z. Let 

-^heat,2 [-^heat,l i ^1 Bpinhorm (3.6) 


with t' ^ |? 7 p ~ (s — 1)^ ^ so we have polar coordinates near ff given by 

l/2(fc-l) 


-p = [{t/P) + + (s - 1)2 + (|y _ y|/x)2) 

^ ■= = 


t X X — X y — y \ ^ ^ 

2^27^ ’ p ’ j ^ith y, z, z. 


x^p 


xp^ 


(3.7) 


Let 

P2 : M2eat,2 ^ M X M X [0, oo)t (3.8) 

denote the blowdown map. Then, similar to the setup at ffi, if we define a = 
(x — x)/x, the interior of ff is a bundle over Y x Ko- whose fiber over p G Y is 
isomorphic to TpY x Z'^ x for T the rescaled time variable below. 

See Remark [Tl] below for further discussion of the need for the two distinct blown 
up faces ff and ffi. 








SPECTRAL AND HODGE THEORY OF ‘WITT’ INCOMPLETE CUSP EDGE SPACES 


15 


(3) Finally, we blow up the time equals zero diagonal, B 2 ■= cl(,d 2 (diag(M°) x {t = 0})), 
where cl denotes the closure, parabolically in time. Note that B 2 intersects the face 
ff at ^ = (1,0, 0, 0), so near the intersection, defining the functions 

_ s — 1 X — X ~ y — y ~ t' t 

V = ^ 


X, 


a = 


Xfc-l 


T = 


x" x'^ 

we have the projective coordinates 


r2{k-l) 


X2/c ’ 


The full heat space is 


x,y,(J,r],T,z,z. 


-^heat “ [-^heat,2j ^2]in/iorrn 


(3.9) 


(3.10) 


(3.11) 


with T ^ 


{z — z)^. The face tf introduced by the final blowup satisfies 
tf° ~ *“T(M), 


(3.12) 


where *'^®T(M) is the incomplete cusp edge tangent bundle defined in (2.8). Con¬ 
cretely, in coordinates (x,y,z) if we set 

yi-yi ^ Za- Za ~k 

Vi = — Ca = — 7r~ 2 : , 




X — X 


(3.13) 


then (x, y, z, ij, C, t) (or (x, y, z, 77, r)) form local coordinates near the intersec¬ 

tion of tf with ff and away from t = 0, and the association ^ dx,Vi ^vi ^ Qa 
induces the map. 

In summary, we have constructed a manifold with corners depicted in Figure]^ 

with a blowdown map /3 as in (1.201, such that has six total faces, three of them being 

the lifts of the standard faces 


If := cl(r'({cc = 0}“)), rf := cl(r'({S = 0}“)), tb := cl(/3-'({t = 0}“)), 


(3.14) 


and then the three faces ffi,ff, and tf constructed (in that order) by radial blow up as 
described above. Denoting the set of the six boundary hypersurfaces by = 

{If, rf, tb, ffi, ff, tf}, and given • G A4(M^gj^j), below we will let p, denote a boundary defining 
function for •, so p, G satisfies that {p, = 0} = • and dp, is non-vanishing on 

• and p, > 0. We can take pg = p as in (3.7) and ps^ = pl~p. Also note that 

a: = PifP = PifPffiPff- (3.15) 

we will show that the heat kernel lifts to be polyhomogeneous on 


In Theorem 
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3.2. Model operators. The blown up space is useful in the construction of a 

parametrix for the heat equation in part because the operator 5* + A (more specifically 
t{dt + A)) behaves nicely at the three introduced boundary hypersurfaces ff,fFi, and tf; in 
particular, the first steps in the parametrix construction involve finding the right asymptotic 
behavior for the heat kernel so that the heat equation (1.14) is satisfied at least to leading 
order at ff,ffi, and tf. 

Thus, we consider the operator A acting on the left spacial factor of M x M x [0, oo)t, and 
the pullback /3*(1(5* -I-A)) to and show that this restricts to an operator at tf. To be 

precise, fix a point p G M and consider the fiber tfp = where tt: tf —diagj\^ = M 

is the projection onto the diagonal (or more concretely it is /3|tf). In the interior of tf, i.e. 
away from the intersection with ff, this is standard [31] . so we concern ourselves only with 
an open neighborhood of the intersection of tf with ff. Indeed, working locally over the base 
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Figure 1. The heat double space (upper left) and the various intermediate 
blow ups together with their blow down maps. 


in both spacial factors, consider a subset of tf of the form {(x, y, z, ry, C, t) : (x, y,T) € O}. 
Now note 

Vtdx = Vidy = dy, = 5^ + 0(t), (3.16) 

and 

^ {rdr - R) , (3.17) 

where R = + rj ■ dy + is the radial vector field on the fiber. 

Letting ttl, tt/j : M x M x [0, oo)t —)■ M denote the projections onto the left and right M 
factors, and End — M x M the endomorphism bundle, whose fiber at {p, q) G M° x M° 
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is End(A*; Ap, for t > 0, the heat kernel restricted to the interior will be a smooth section 
of this bundle. To study the heat kernel at the boundary we use the incomplete cusp edge 
forms and the corresponding endomorphism bundle End(“^®A*) back to M x M x [0,oo)t 
and then to via the blowdown /?. As usual, restricting to the spacial diagonal gives 

the ‘little endomorphism’ bundle 

EndrA*)Uiag(M) endrA*) 

where end(“’®A*) —M is the endomorphism bundle of the exterior algebra of M. The 
restriction to the time face, (3* End|tf, is isomorphic to the pullback of end(Ap to the tangent 
bundle of M via the projection map. 

Writing w = {x,y, z),w = {x, y,T), sections of 13* End near the fiber of tf over p can be 


written 

a = aijdwi 0 dwj, (3.18) 

i.J 

where I, J run over all multi-indices and dwj is dual to dwj, and here ajj = aij{w,w,t). 
We claim that, writing sections of (3* End near tf as sections of (3* Endjti — A*{M) 0 A(M), 

t{dt -I- A) = (i {rdr - R) + ct(A)) 0 Id-I-O(T), (3.19) 

where (t(A) is a constant coefficient differential operator in the coordinates S = 
depending on the metric g at p = {x,y,z), namely 

a{A) = {ds+Xg^p^dB.Xg(p))‘^, (3.20) 


acting on differential forms on the vector space '“ApM) with metric g{p). Indeed, let w 
be geodesic normal coordinates. In the interior of tf away from ff we have coordinates 
E= {w — w)/^/l, w, y/i. Then t{dt + A)a = {t{dt + A) j aijdwi) ® dwj and moreover 

-kdwj ® dwj = {-kdwj) ® dwj 

= zLdwin 0 dwj -b 0{w — w) 

(3.21) 

= ±{dwia + y/idEja) 0 dwj + 0{ViE) 

= {*g(p)dwi) 0 dwj + 0{t), 

Similarly, letting the exterior derivative act on the left gives d(adwi 0 dwj) = {d^^adwi A 
wj) 0 dwj. 

To motivate our construction of the heat kernel further, in a neighborhood of tf let 7 be 
a section of End with the property that 7 I diagj;^ = Id on the form bundles, and consider 
the section of (3* End on of the form 

such that G{p,q) satisfies that j3*{G{p,q)/t)\ti = ||2||^, that is, that G{p,q)/t restricts to 
the metric function on tf. Such a form 7 and function G can be constructed but we neither 
prove nor use this; we merely use it as motivation. It is straightforward to check that for 
any smooth compactly supported form a 

lim f K{p,q,t)a{q)dYo\q = j e“'^“'' 3 <f’)^^a(p)\/( 7 (p)|dS| = a{p), (3.23) 

Jm jm 

and in fact the convergence takes place in L^. (In fact, such an endomorphism 7 can be 
constructed easily by taking the identity map on over M, pulling this back via (3 to 
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/3* Endjtf and extending off smoothly in a neighborhood. On each exterior algebra A*M, 
the identity can be expressed in terms of a basis with dual basis e* as ei ® e*j. In a 
neighborhood of tf n ff we can take the basis Ci to be dx, dj/i, x^dza, i.e. we can take the Ci 
to be a basis of forms for '“A* all the way down to x = 0.) 

Working in the coordinates (3.13), since the Taylor expansion of the 

heat kernel at tf should take the form 


1 




j=o 


(3.24) 


where the bj = bj (a?, y, z, ^,r],C) are sections of (3* End, which we again write in a neighbor¬ 
hood of tf n ff as sections of End(*“A*) pulled back to the fibers of tf. Writing each bj as a 
finite sum of terms of the form 

a®g~^j3, (3.25) 

where a and f3 are sections of *'=®A* and g~^ indicates taking the dual vector field, we see 
that by (3.19) we have. 


1 


(a - ® Id &o = (a - a-R + 


1 


Ah 

0 


0 


) (g) Id ) 6o 


(3.26) 


The only solution to this equation which gives the identity operator at t = 0 is 

= g-||H||V4 X Id. (3.27) 

The procedure of solving for the remaining bj is standard m Chapter 7]; letting the 
Laplacian act on this expansion we show that on each term Uj it acts fiberwise like a 
constant coefficient, second order elliptic differential operator plus the radial vector field 
plus a constant corresponding to the order of the term in the expansion. We have the 
following 

Lemma 8. There exist sections bj o/yiphg(Endjtf) satisfying 

bj = e"ll“ll'/%-(x,y,z,^,r?,C) 

where bj is a polynomial in g, f and a polyhomogeneous section of End over tf, such that 
for any distribution H' in ^P*'®(End) with asymptotic expansion near tf given by (3.24) we 
have 

t{dt + A)i7' = 0(r~), 

i.e. t{dt + A)H' vanishes to infinite order at the blown up t = 0 diagonal, and, moreover, 
the asymptotic sum of the bj exists and yields such an H'. 

The existence of a distribution H' as in Lemma is only a first step in constructing a 
parametrix for the heat kernel. We will discuss the rest of the process in §3.3[ 

A useful double check of the order of blow up of the heat kernel at ff is the following. 
Near ff n tf we have 

5{x — x)S{z — T)S{y — y)= 5{fiTx^)5{gTx^)6{C,T) 

Since Id = lirntx^o H dVolg ^ lirntN^o Hx^^dxdyd'z, we confirm that H should have order 
—nk at ff. In fact, we can deduce more; considering x^‘^H\ff, on the interior of ff we can 
use coordinates in (3.10), we get that 

S{x — x)5{y — y)S{z — z)= 5{a)5{rf)5{z — z), 


(3.28) 
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which means that, on the face ff, we expect that the restriction will be given by 

S{a)6{rj)S{z — z) at least as the time variable T = t/x^ goes to zero, as that is the region in 
which the action of H is definitively approximated by the identity. On the other hand, x 


commutes with the heat operator dt + A. As we will see in (3.31), t{dt + A) restricts to an 


operator on ff and defines a fiber-wise heat type operator on ff, so we expect to have 

t{dt + AMx^>^H)\e = 0. (3.29) 


This, together with (3.28), implies that an ansatz for the heat kernel should include that on 
each fiber of ff, x^^H[g is the fundamental solution to the induced heat equation 
on the fiber, more precisely, it is the solution which equals Scr=o^rj=o Idz at time equals 
zero. The induced heat equations are translation invariant in a and rj, thus induced by 
convolution operators, and the heat kernels we speak of are the convolution kernels in cr and 
V- 

As for the blow up at ffi, as we will see below, the operator acts as a modified heat 
operator in and Y on the bundle of fiber harmonic forms, so in the coordinates in (3.5) 
we will have 

5{x - x)5{y - ^S{z -z) = - l)<5(z - (3.30) 

In this case, t{dt + A) only admits a restriction to ffi on the fiber-harmonic forms H, on 
which S{z — T) becomes projections nz,y onto the kernel of A^’^. Thus we expect that 
on fiber harmonic forms is given by the convolution kernel for the heat kernel 
in ry, times the dilation invariant kernel for the heat kernel in s with limit <5s=i at time 0. 

We now compute the asymptotic behavior of t(dt + A) at the faces ff and ffi. First we 
will work at ff. 


Proposition 9 (The model problem on ff). The operator Ng{t{dt + A®)) = t{dt 
acts fiberwise on ff, and is expressed in the coordinates in ( M by 


Ng{t{dt + Aa))=T{df + 


-dl + A,, -f A^’« 

0 


0 

-dl + A„ + A^^y 


A9)|f 


(3.31) 


on the fiber above y € Y. Here A^ is the constant coefficient Hodge-Laplacian on the 
tangent space TyY with translation invariant metric h{y), and A^’^ is the Hodge-Laplacian 
on (Z, ky) 


The situation is more delicate at ffi. As we will see in Section [3i3l near ffi, it will suffice to 
consider t{dt-\-A) restricted to approximately fiber harmonic forms, (see Section 2.3). Thus 
let 7 e ai®'Hphg(‘“A*) with the space on the right defined in (2.28), and thus by assumption 


S^'y,d^'y are both 0{x^~^'^). From (2.15) it follows that for such fiber harmonic forms 

Ay = Ao7 + x-^’^Azj + x-\6zQi + dzQ 2 )l + 0{x^-^) (3.32) 


-2kx-'^-^d^ 
k{f -N){1- kN)x-^ 


where Aq acts on forms decomposed as in ( |2.6| ) as 

kf. 


Aq — — dj, - dx + Ay -f 


fcN(l-fc(/-N))a;-2 

-2kx-^-^5^ 


By fiber harmonicity, x~'^^Az^ = 0{x^~^). Thus the two terms x~^^Az and x~’^6zQi act 
on fiber harmonic forms as operators of order x~^, and thus in the heat operator t{dt -\- A) 
there are term behaving like tx~^ (on fiber harmonic forms) and t/x~^ is not a bounded 
function at ffi! On the other hand, if we project back to the fiber harmonic forms, then by 
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(2.31)-(2.32) we kill these terms; concretely, with II-h the fiber harmonic projector in (2.30), 
we have 

n«An« = Ao + x-^E' (3.33) 

where E' G Diff^ pj^g (see (2.13l), and thus does not decrease the order of vanishing of 
polyhomogeneous distributions. Defining 


Pa,b ■— - ds H—5- 


and 


(3.34) 


(3.35) 


a(N):=fc/, /3(N) :=fcN(l-fc(/-N)), 7(N) = fc(/- N)(l - fcN), 

we have the following. 

Proposition 10 (Heat operator on fiber harmonic forms at ffi). Restricted to the fiber 
harmonic forms H as defined through (2.25), 

NeMdt + ^^)) :=n«t(at + A)n«|ff, (3.36) 


restricts to the face ffi in the coordinates ( |3.5[ ) as 

NsAt{dt + Aa)) = t' 


dt> + 


a(N),/3(N) 

0 


0 




■ A„ 


(3.37) 


a(N).7(N) 

Remark 11 . Analysis of the fiber harmonic forms is necessary in particular because the 
structure of the operator A® is such that, off of the fiber harmonic forms, the leading order 
term is x~‘^^, while restricted to the fiber harmonic foms the leading order term drops 
in order. Indeed, if it weren’t for the presence of the term x~^P in (2.15), which presents 
complications in the analysis, on fiber harmonic forms A® would be given by to leading 
order by Ag. Indeed, the need for the two different regimes represented by the boundary 
hypersurfaces ff and ffi is exactly this change in asymptotic order of the operator on and 
ojf the fiber harmonic forms. Correspondingly, we will see below in the proof of Lemma [7^ 
that the operator t{dt + A) restricted to ff has a fundamental solution which vanishes at ffi 
to infinite order off the fiber harmonic forms. 


The heat equation for the regular singular ODEs in (3.34) has been studied in detail. To 
such an operator there corresponds a pair of indicial roots given by the order of vanishing 
of homogeneous solutions, specifically Pa,b{s^) = 0 if and only if 


1 = 


-{A - 1) ± ^(H- 1)2 +4H 


(3.38) 


The numbers £ give important information about the operator Pa,b, in particular they give 
the order of vanishing of the Green’s function at s = 0. The operators that will arise in our 
work are those in the matrices in (3.37). We define the indicial set 

/ 

A=U 


— {a — 1) ± ^/(a — 1)2 + 4/3 —(a — 1) ± 1 /( 0 ; — 1)2 + 47 


N=1 

/ 


(3.39) 


= U - l)/2 ± |fc(N - //2) + 1/2|, -{kf - l)/2 ± |fc(N - //2) - 1/2|} 


N=1 


Letting 


=B + 


,A-1 


)^ > 0 


(3.40) 
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where u > 0, from [35l Vol. 2, Eqn. 8.45] there is a fundamental solution Ha 

{dt + PA,B)HA.B{s,s,t) = 0, and H ^ Id as t —)■ 0, on L^{s^ds). (3-41) 


Indeed, one has the explicit formula 




ss 


2t 


(3.42) 


where is the modified Bessel function of order u of the first kind [1] Chap. 9]. 

As discussed below (3.30), at the face ffi we expect the heat kernel to be of order 
~-i-b-kf ^ Thus we expect to have 

1 


0 = tidt + A)H=^^^^{t{dt 


.A))(ji+b+fe/i/), 


(3.43) 


and since n^t(9t + A)n>^ defines a differential operator on section oi H ^ H restricted 
to ffi, w e in clude in our ansatz for the fundamental solution (1.14), and indeed prove in 
Theorem 


12 


below, that there is a fundamental solution H satisfying that H has a 

smooth restriction to ffi, and writing 

NedH) := H)\s„ we have NsMdt + A))NsAH) = 0. (3.44) 


Furthermore, again as discussed below (3.30), it is sensible to include in the ansatz for H 
that Ng^{H) is the fundamental solution for the model operator Ng^{t{dt + A)), meaning 
specifically that Ng^ {H) is a section of the restriction of the sub-bundle End('H) to ffi and is 
given using the fundamental solutions to the model heat equations Ha,b from ( |3.41[ )-(3.42 ). 
Specifically, we will have as an ansatz that Ng-^^ (H) = Kg ^, where 


■ — 


/ i7„./3(s,l,C) 


0 


l,f') 


(47rt') 




(3.45) 


where a, [3,^ are as in (3.35), and in particular continue to be operators depending on the 
fiber form degree N. The distribution Ng^{H) is polyhomogeneous on ffi, and the leading 
order behavior at s = 0 satisfies that for 0 < c < < C < oo, for some smooth a{t'), b{t'), 

1, t') - s-('=-^-i)/2a(t')s‘'^“’'^\ i^a.7(s, 1, t') - (3.46) 


with V as in (3.40) 


r fc(//2-N)-l/2 ifN<//2, 

\ fc(N - //2) + 1/2 if N > //2, 

/ fc(//2 - N) + 1/2 if N < //2, 
\ fc(N-//2)-l/2 ifN>//2, 


and thus by (3.46) on ffi in the region 0<c<t'<C< oo, 

-kN 

Kffj = 0{s'') where T(N) = —kN + 1 

-fc(/-N) 


if N < //2, 
if N = //2, 
if N > //2. 


(3.47) 


(3.48) 


In words, each has two indicial roots, the order of Ha,p for fixed s, f > 0 is the larger 
of these two, and p is the smaller of the leading orders of Ha^p and 77ct,7- 

The behavior of the heat kernel at ffi also shows what to expect at the left face, the lift 
of X = 0. There we should just have the projection onto the fiber harmonic forms times the 
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®A*, times the lifted 


~ K := 


Hy 


(3.49) 


leading order behavior of the Ha^p and acting appropriately on 

heat kernel of the base Y. Indeed, we expect 

HaCM),pCN){x,X,t) 0 

0 '^a(N),7{N) 

where Hy is the heat kernel on (Y, h) lifted to the tubular neighborhood U in (2.1) via the 
projection and k acts on sections of the bundle of fiber harmonic forms % with its grading 
by hber form degree N. (See Section 2.3 ) In fact, with V = i?(N) the fiber degree dependent 
weight in (3.48), 

K e X M;04=o'^’^ (3-50) 

where % and %, respectively, the pullbacks of the fiber harmonic form bundle (defined on 
a neighborhood U of the boundary) via the left and right projections of M x M. 

As discussed below (3.29), on the face ff, we expect that the heat kernel will have leading 


asymptotic x 


—nk 


so we expect and prove that 

Ns{H) := (2?”'=if)|ff, ^ Nsitidt 


A))NsiH) = 0. 


(3.51) 


Again, we will set Ng{H) equal to a fundamental solution to the heat equation, namely, 
using the decomposition in (3.37), we expect to have Ng{H) = Kg where 

Kgja, 77, z, z', f) := Id2x2(47rf (3.52) 

where Hz,y = Hz,y{z, z',T) is the heat kernel for A^'^. 


Before stating the full structure theorem for the heat kernel let us briefly recall the notion of 
an index set, which by definition is a set of exponents £{•) = {( 7 ,_p)} C C x N associated 
with each face • G {If, rf,tb,ffi,ff,tf} such that 

(i) each half-plane Re ■j < C contains only finitely many 7 ; 

(ii) for each 7 , there is a number P{'f) G No such that ( 7 , 73 ) G f (•) for every 0 < p < 
Pil) and ( 7 ,p) ^ £(•) if p > P{-f); 

(iii) If ( 7 ,p) G £{•), then (7 + j,p) G £{•) for all j G N. 

We give a full rigorous definition of polyhomogeneity in Section but roughly speaking, 
we call a differential form a polyhomogeneous with index family £ = {£(•)!•€ 
{If, rf, tb, ffi, ff, tf}} if it has an expansion at each boundary hypersurface • with exponents 
determined by the corresponding index set £’(•) and coefficient functions which are them¬ 
selves polyhomogeneous (with exponents determined hy £). For example, smooth functions 
on are polyhomogeneous with indicial set satisfying £{•) = Z x {0} for all •, and if a 

polyhomogeneous function vanishes to infinite order at a particular boundary hypersurface 
•, then it is polyhomogeneous with an index set £ satisfying £{•) = 0. We define 


inf£:(*) = inf {7 | {y,p) G £(•)}. 


Theorem 12. There exists a section H of /?* End over which is polyhomogeneous, 

i.e. H G x4pi,g(M^ggj; /?* End) for some index set £ and satisfying the following properties. 

(1) In the interior of {dt -f A)H = 0. 

(2) The operator Ht defined initially on forms a G CIP(M;fl*(M)) by 

Hta{w) = / H(w,w,t)a{w) (TYo\^ (3.53) 

Jm 

is symmetric on L^(dVolg), and for such a 


Hta -G a as t ^ 0 


( 3 . 54 ) 
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in . 

(3) The index set S of H satisfies that f (tb) = 0, while 

inf £(ffi) > -1 - 6 - kf, inf £(fr) > -kn, £(tf) = N - dim(M), (3.55) 

where £{•)>€ means that inf^ pg 5 (,) Re C, > c. 

(4) Moreover, at the faces ff and ffi, 

{x^+^+>^f H)\s, = = Kff, (3.56) 

with and Kg the model heat kernels defined in (3.45) and (3.52). 

(5) At the left face If, with k as in (3.49) 

= /c(l + 0(pff)). (3.57) 

In particular, for t > 0 fixed, Ht is approximately fiber harmonic in the sense of 
(2.25), and p\f is a boundary defining function for li (see (3.15[)J. Thus 


inf £ (If) > —Y + 1- 


(3.58) 


The behavior of Ht at the right face rf can be deduced from symmetry. Moreover 
for the behavior at the codimension 2 face If H rf, the leading order behavior is the 
product of that at If and rf, i.e. Ht = 0((pif/3rf)~^^^^''’^)- 


The proof of Theorem is at the end of ! 

3.3. Parametrix construction. We will establish the following 

Proposition 13. There exists a poly homogeneous section 

if e.Aphg(M2eat;/3*(End)) (3.59) 

such that K satisfies (3.54)-(3.57) above, and such that Q := t{dt + A)if, which is polyho- 
mogeneous, has index set £' satisfying 

£'{Si)>-l-b-kf + l, £'{S)>-kn + l, £'(lf) = £'(tf) = £'(tb) = 0. (3.60) 

Our work above allows us to break the proof of the proposition into two main steps; first 
we will prove the following lemma: 


Lemma 14. There exists a poly homogeneous Ki satisfying (3.54)-(3.57) above, together 
with (3.24) for the indicated bj. Furthermore, Ki can be taken fiber harmonic in a neigh¬ 
borhood of ffi. In fact we can assume that ifi = k in a neighborhood of If and also in a 
neighborhood o/rf, where k is as in (3.57). 

Proof of Proposition \l 3\ assuming Lemma \14\ Assuming that we have such a distribution 
ifi, we study t{dt + A)ifi. Automatically we have that t{dt + A)ifi vanishes to infinite 
order at tf and tb, as follows from Lemma Furthermore, t{dt + A)ifi vanishes to order 
—kn + 1 at ff by (3.31) and the fact that the leading order term Kg there solves the model 
problem. 

At ffi things are again more delicate. Recall that ifi = at ffi, where is 

the boundary defining function for ffi in, e.g. pg^ = pl~p with p as in (3.3) and p as in (3.7). 
Since ifi is fiber harmonic near ffi, by (2.30) and (3.32[) we have 


dzQ2)Ki 


Aifi = Aoifi + x-^'^A^’^ifi + x-^{dzQi 

= Apifi + X ^{6zQ\ + dzQ2)Ki 
+ x-\dzQi + 5zQ2)K' + 


x-^EKi 


■x-^EKi 
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Furthermore, by (3.33) we have that + ^)^uKi is order —1 — 5 — fc/ + l since its 

leading order term solves the model problem. 

We assert the existence of a polyhomogeneous distribution A of order —1 — h — kf + 
k such that t{dt + A)(iFi — A) itself vanishes to order —1 — b — kf + 1 at ffi. Indeed, 
since the leading order term in t{dt + A) is and since by (2.33) we can solve 

/\Z,yB = (dzQi + SzQ 2 )K' + SzQiKi + where B is polyhomogeneous with 

asymptotic expansion determined by the expansion of the right hand side, in particular 


B = 0{pg^ ^ We take A = x^B and thus obtain, with P as in (2.15), 


t{dt + Ai)(Ki — x^B) = t{dt + Ao)(Ali — x^B) + tx ^E{Ki — x^B) 

- tx-’^Px^B + 

= t{dt + Ao){Ki - x^B) + tO{pQl~'"~^^) 


iO{Ps^ 


+ tx 

Since the expansion of B at ff has the same order as Ki, the distribution 

K 2 = Ki- x^B 

has all of the desired properties of K in the statement of the proposition except that (9* + 
A)K 2 is not rapidly decreasing at If. Note that, since Ki = where V is the 

(fiber degree dependent) order of n computed in (3.48), by well-posedness B also satisfies 
B = 0{pl). 

To deal with the expansion at If we argue along similar lines, but there we iterate the 
argument to get a parametrix K with {dt + A)K vanishes to infinite order at If. (We work 
in the interior of If though the arguments at the intersection of If and ffi are the same in the 
projective coordinates s' = xjx^rl = {y — ’dpjx^T' = tjx^ together with z,x,y,'z.) Recall 


that Ki = K near If and thus K 2 = k — x^B near ffi. Again with P as in (2.15), we have 


{dt + A)K2 = x-^Pk + x'^-'^Ek - x-^A^^yp - x-’^Px'^B + 0{x''+^) 
= x-^{dzQi + 6 zQ 2 )k - x-^A^’yp + 0(x"+'=-2), 


(3.61) 


where v is the leading order power of k computed in (3.48). As in the argument at ffi, since 

there 


2.3 


the RHS of (3.61) manifestly gives that YVuPdt + A)K 2 ) = 0{x''^^~‘^), by Section : 
is distribution Aq such that x'^A^’^Aq = {dzQi + 6 zQ 2 )k — A^’^R + 0(a;'^+^). Here the x'^ 
in front makes it so that Aq is 0(1). Thus 

{dt + A){K 2 - x^+^AP, = 0{x^+'^-'^) - x-^Px^+'^Aq = 0{x^). 

We will now solve away iteratively to decrease the order of the error. For this we assume 
for the moment that we are given, for some q > V + e, any distribution Ai = x^Ai + 0 {x'^'^’') 
with Ai smooth and non-vanishing up to the boundary as in ice—form. First, we find a 
distribution Ri so that a;'^A 2 := {dt + A)(a;'J+^^Ri) — Ai is fiber harmonic. We can do this 


by solving (/ —n^)Ai = A^’^Ri -I-O(x^) as in Section 2.3 where H^ is the projection onto 


the fiber harmonic forms, since then {dt + Ajx^+^^Ri = x'^A'^’^Ri -|- 0(x^+^). We then 
construct a term Ci with {dt + Ajx'^+^C'i « A 2 , as follows. Decomposing A 2 = (A 2 , A|) 


as in (2.6) in the left varibles, and noting that if Ci = {{—{q + 2)^ — {a — l)(g -I- 2) -|- 
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> (-(? + 2 )^ -{a- l){q + 2 ) + 7 ) ^A^), then 


Pa,0 

0 


0 

Pn ' 




(The numbers we divided by above are non-zero, since the indicial roots of Pa,p and -Pa ,7 
are bounded above by T — e, as explained below (3.48).) For this Ci we have 

x^A 2 - {dt -f A)a;«+2C'i = - Aox'J+^C'i -h x-’^P'x'^+'^Ci + 0{x‘^+^) 

= 0(a;«+'^) -f x-'^P'xi+^Ci + 0(a:«+2+'=-2), 

S can be taken to be the order of the subsequent term in the expansion of A 2 

we see that the left hand 


where q 

where Aq is in (|3.32|) and P is as in (|2.15|), and thus by Section 


2.3 


side lies in the image of to order x'^. We can thus find a distribution Di such that 
x‘>A 2 - {dt + A)(a;«+2Ci - cc'^+^+'^Pi) = 0{x‘^+^) - x‘>+‘^-^A^^y+ x'^P^x^+'^Ci 

= 0(a:«+i), 

which gives 

{dt A)(x«(a;2'=Pi - x^Ci + x^+'^Pi)) = x^A^ + 0{x’^+^). (3.62) 

It is straightforward to check that the added terms do not increase the order of blowup at 
ffi. Thus we can kill off the leading order term of x^A, and in fact can kill off all terms 
iteratively by this process. (If there are log terms present the arguement is analogous and 
left to the reader.) 

From the previous two paragraphs, it follows that we can find a distribution K' such that 
K := K 2 — K' satisfies the requirements of the lemma, specifically such that t{dt + A)K, in 
addition to having the same leading order asymptotics at tf and ff and ffi that t{dt + A)K 2 
has, also vanishes to infinite order at If. Indeed, since we can solve away terms to obtain 
errors of succesively decreasing order, taking the Borel sum [29] of these distributions gives 
K'. □ 


Proof of Lemma If We seek a distribution K' with the stated asymptotic properties; namely 
(3.54)-(3.56) and (3.24) for the indicated bj. Such a K' will exist by Lemma[2l]in Appendix 


|A| provided the hypotheses are satisfied, meaning that the following matching conditions 
hold. We must find a set {p,} of boundary defining functions for the boundary hypersur¬ 
faces, • = If, rf, tf, tb, ff, ffi of such that 


Kff = 


(47r)"/2 






(3.63) 


Kfi on ff nffi. 


and that ks, Kg-^ and the bj vanish to infinite order at tb. Indeed, in the notation of Lemma 
we have ki = {pgiKg^ and K 2 = {pg/x)^'^Kg, and the matching conditions in 


21 


terms of ki and K 2 in Lemma 21 are exactly (3.63). We use boundary defining functions 


Ps = PiPSi = pIp for the faces ff and ffi defined in (3.7) and (3.3). We define boundary 
defining functions pif,Prf for the faces If and rf by the equations 

PifPffPffi = X, prfPffPffi = X. (3.64) 


Finally, we use r in (3.13) as ptf, though it is not valid at tbntf, all the distributions in 


question will vanish to infinite order there and there will be no conditions to check. 


The first matching condition in (3.63) follows easily since the coefficients of the expansion 
of are determined by the same differential equation which determines the bj, and the 
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coefficients in both expansions are uniquely determined by their being equal to polynomials 
times Gaussians on the fibers of tf n ff. 


Finally we check that the second condition in (3.63) holds. First we consider = 
above the point y €Y (i.e. restricted t o fFi„ ). In the polar coordinates in 
(3.6) and using the boundary defining functions above (3.64), we have 


s = -h 1, 


t' = 


Pli 


V = 




k -1 


Pit Pli Pit 

Using [U Eqn. 9.7.1], we have that the modified Bessel function satisfies 

/y(z) = (e“"'/V27rz)(l + 0(1/2:)), 

and thus 


(3.65) 


PsiPitPft) 


-i-b-kf 




-kf 

Pit 


(47r,/jU+i)/2 


e-(IV'.+lbd5)/4'/'t(i + o(pff)). 


(3.66) 


On the other hand, above each base point y gY, K.gy{a,ri',z,z',T) can be written using 
separation of variables with respect to the spectrum of A^’^. Indeed, since Hz,y has discrete 
spectrum, it is standard that Hz^y{z,z,t) = Ho + i? where Bq is projection onto the kernel 
of and \E\ < e~^°* as t —>■ oo, Aq being the smallest non-zero eigenvalue of Az,y Thus 


Kff^y = ( 27 rT)-('’+l)/ 2 g-(<T"-FVI^,)/ 2 n^ 


where E' is exponentially decaying. Now we have 

V' = 


rji _ 1 2 (fc 1 ) — 2 k 

T — (PtPsi Pit ■> 


-ik-l)-k 

Pit y 


and thus 


Pst^^'^HpitPfti 


,—kn . 


Kft = 


-kf 

Pit 


T ~(k — l) —k 

t^ = V’xPs, Pit y 




(3.67) 


(3.68) 


(3.69) 


(47r(/jU+i)/2 

so the matching condition at ff n ffi holds. 

For the behavior of Ki at If, since Kg^ is given by n at the boundary and we 

can take Ki to be equal to k, at If and at rf. □ 

We can now conclude the proof of Theorem [T^ modulo arguments in Appendix 


Proof of Theorem 12. Having established the existence of a parametrix K as in Proposition 


13 we will now prove Theorem 12 We do so by inverting the error Q = t{dt + A)K 
from Proposition |13| via a Neumann series. To be precise, it will be convenient to think of 
distributional kernels A{p,p',t) on M x M x M+ acting on C^{M° x (0, oo)) by operating 
as convolution kernels in the time variable, so for (j) G Cf°{M° x (0,oo)) by 

{A-k (j)){p,t) '■= [ f A{p,p',t — s)(j>{p', s)ds dVolpi. (3.70) 

Jm Jo 


Then 


{dt + A)K = I + t-^Q 


(3.71) 


and the right hand side can be inverted via a Neumann series, i.e. (Id +t ^Q){I + Q') = Id 
where Q' = ^“^d {t~^Qy = t~^Qk ■ ■ ■ kt~^Q, j-times. In Proposition 


28 below, we show that the summands {t ^Q)-' vanishes at successively faster orders at ff 
and ffi. Moreover, as discussed in [23 E], this series is convergent in C°°, and the infinite 
order of vanishing of t~^Q at If is preserved in the sum, i.e. Q' vanishes also to infinite 
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order there. In fact, one sees as in m that K{I + Q') is polyhomogeneous with the index 
set £ statisfying the properties of Theorem □ 


4. Spectral and Hodge theoretic properties of the Hodge-Laplacian 

In this section we deduce the main theorems from the introduction. We begin with a 
detailed analysis of the polyhomogeneous forms in the maximal domain. 


4.1. Polyhomogeneous forms in I?max and Recall the definition of 'Dmax and 

'Dmin from the introduction, and the space .Aphg(''^®A*) of polyhomogeneous ice-forms (below 
denoted simply by -Aphg) discussed in Section]^ We determine conditions which assure that 
a given polyhomogeneous differential form 7 € -Aphg is contained in the maximal domain 
I^max of A®. This will be used to show, with an additional assumption on the index set of 
a phg form, that 

T € ^max 14 -Aphg 4 7 ^ ^^min 14 -Aphg. (^R) 

Let 7 G ^phg be contained in the maximal domain, i.e. we assume that 7 G and 
A ®7 G L^. Let 7 = where 7 = 70 ( 2 /, 2 ) + 0(x'^). Here notation such as 0(x'^) indicates 
that the differential form 7 is locally a combination of basis forms 

dyi A and dx A x^^dyi A x’^ ^dzA, 


where I and A are multi-indices on the base and fiber, respectively, with coefficient functions 
which are bounded by cx'^ pointwise in norm when a; 0, and 70 is a form on M whose 
coefficient functions are independent of x. Let us determine the possible range of values 
s. From (2.9) it follows that in a neighborhood of the boundary, the volume form of the 
cuspedge metric g is 

dVolg = x’^^pdx Ady A dz, 

where p = a{y, z) -\- 0{x^) and a is a non-vanishing positive function. It follows that 


x^j G L‘^{M,g) 


s > -^(fc/ + l)- 


(4.2) 


We begin by analyzing the indicial roots of A®, specifically we find the order of vanishing 
of approximately fiber harmonic homogeneous forms in the kernel of A®. By Proposition]^ 
the leading order part of A® restricted to approximately fiber harmonic forms (see Section 


2.3) is 


n«A®n« -:= 




a(N).,3(N) 

0 


0 




a(N).7(N) 


with Pq(n). 7 (n) the operators, depending on fiber degree, defined in (3.34)-(3.35) 

We note that 


-Pa(/-N)./3(/-N) - -Pa(N).7(N) 


(N = 0,...,/). 


(4.3) 


By (3.39), the values s for which Pq(n),;S(n)®^* = 0 for which elements of size x^ are in 
addition contained in L^{M,g) is 


-fcN 

I-fc(/-N) 


if N<i(/+i), 

if N>i(/-i). 


(4.4) 


Proposition 15. Suppose the differential form 7 = ( 7 ^, 7 ^) = G Alphg and 

that 7 ® = fQ{y,z) + 0(x‘^) is contained in the maximal domain I?max- (Thus the leading 
order term is assumed not to have a logarithm, as is a priori allowed for phg-distributions.) 
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Then each Sj is an indicial root o /for some 0 < Nj < f or Sj > li-kf + 3). 
In either case, Sj > ^(—kf + 3). 

Proof. Recall form Proposition that 


AS = Ao + x-^P + x-^E 


where 


An = 


P 


a(N),,S(N) 

0 


0 


p 


q(N),7(N) 


^A^’S + Ay 

—2kx~^~^5^ 


-2kx-’^-^d^ 

1 + Av 


(4.5) 


In view of the symmetry (4.31 it suffices to consider the image of the the component 7 ^ = 
2 ;Si^i under A®. The discussion naturally falls into several cases. 


(1) The form 7 q is not approximately fiber harmonic in the sense of (2.3 Then the 

lowest nonvanishing term in (4.5) is A'^yg, which is contained in if and 

only if 

Si > - !)• 

(2) The form 7 ,) is approximately fiber harmonic harmonic. We then consider the 
following subcases. 


(2.a) Si is an indicial root of Pa(Ni),/ 3 (Ni) and hence equals the number in (4.4). 
(2.b) Si is not an indicial root of Pa(Ni),/ 3 (Ni), he. Po!(Ni),, 9 (Ni)(j^*’^ 7 ^) 7 ^ 0. We 
claim that at least one of the following two statements holds true: 

* The polyhomogeneous expansion of 7 ^ contains a term 7 ^ of order 0{x^) 
where 5 — 2k < si — 2 and 7 ^ is not approximately fibre harmonic. 

* The lowest nonvanishing term in the first component of A ®7 is of order 


If this claim holds true we conclude that the lowest nonvanishing term in the 
first component of A ®7 is of order at most To prove the claim, assume 

that the first statement is false. Then the second one must hold true as is clear 


from the form of the Laplacian Ag in (4.5). To be specific, collecting the terms 
of order in the first component of Agy we obtain 




-'^dzQ2T^ + x-’^QidzT^ 


(4.6) 


f’a(Ni)./3(Ni)(a;'''7^)+ a; ^'“A^ 

+ x-HzQit^ Px-'^QsSzt^ 
for suitable differential forms r^,..., r® of orders 
ri = 0(x®i+2'=-2), = Cl(x"i+'=-^), and r^' = 0(x"i+'=-2) (j = 3,4, 5, 6 ). 

By Hodge theory, the terms dzx"^ and Szt^ both vanish approximately in the 


sense of (2.3 since otherwise a nonvanishing term x respectively 


X would occur. These are both of order strictly less than si — 2, 

contradicting our initial assumption. Considering the remaining five terms in 


(4.6) it follows from Hodge theory and the assumption that 7 q is approximately 
fibre harmonic that the sum 




d^T^ + X-^dzQ2T^ + X-^SzQiT^ 


(4.7) 
Hence 

we conclude that the nonzero term Pa(Ni),,s(Ni)(3^'’^7^) cannot cancel with the 
sum (4.7). It follows that the second statement is true, whence the claim. 


is approximately orthogonal to Pq;(Ni),/ 3 (Ni)( 3 ;'*’^ 7 ^) in the sense of (2.3 
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The asserted statement follows by inspection of each of the above cases. In case (1) it follows 
from 

Si > ^(3^/ - 1 ) > + 3)) 

using that k > 3. In case (2.b) the lowest nonvanishing term in A ®7 is of order at most 
Si — 2. Since 7 S T’max it follows from (4.2) that 


Si — 2 > —-{kf + 1) 


Si > —-{kf + 3). 


In case (2.a), the form 7 q is approximately fibre harmonic and therefore by the Witt con¬ 
dition N ^ The exponent si is given by (4.4) from which it follows that if / is even 
that 

Si >-fc({-l) 

(here we use the Witt condition) and if / is odd that 

Si > -k{- - ^) > - — 

12 2 ' ~ 2 


3 

2 ’ 


where the last inequality follows from the assumption k > 3. 


□ 


Lemma 16. Assume k >3. Then Umin H .4phg = T’max H .Aphg. 

Proof. It suffices to prove the inclusion Umax H .Aphg C I?min H .4phg. For e > 0 we define 
the logarithmic cutoff function : [ 0 , 00 ) —)• [ 0 , 1 ] by 


0 , 

Y (r) ■- J log(a;/e^) 

XeWi •- ^ log(e) I 

1 , 


For < X < c it satisfies 

X(.(x) = and 

log(e)x 

Let 7 G Umax n yiphg and set "/e = Xe"f- Then 


X < e^, 

< X < e, 
X > e. 


log(e)x 2 - 


(4.8) 


A®7, = x,A97 - {dlxeh - (5.Xe)(5.7) - -fid.Xeh 

+ {d^Xe){dz,x) + B\da;Xe){dy,x), (4-9) 

where A^ = O(x^) and B’’ = 0(x^^) are bounded functions with that order of decay in x. 
We show that 


||A® 7 e - A®7||i2(M-,g) 0 

hence establishing that 7 G l?min- It is clear that 

WXe^^l - ^^l\\L^M,g) ^ 0 


as 


0 , 


(4.10) 


as 


0 , 


and thus it suffices to consider the next three terms in (4.9) and to show that 
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converges to 0 as e —>■ 0. Let 7 = a ;®7 for some 7 = 0{1). A short calculation shows that 
each integrand in (4.11) is of order for some i5 > 0 and hence converges to 0 as e —>■ 0 

if 

kf , 3 


s > —- 


(4.12) 


In the borderline case s = — ^ + | we still get convergence since then the first integral in 
(4.11) becomes 

1 


log (e) 


— dx = 


^ (log(e) - log(e 2 )) = - ^ 


log (e) 


log{£) 


0 as £ —)■ 0 , 


and analogously for the second and third integral. Hence 

^ 7 e I?n 


^^-T + 2 


for any 7 = x ®7 G I?max H Aphg- On the other hand, Proposition 15 shows that 

7 = X 7 ^ 2 ?niax O ^phg — 




and hence the claim follows. 


□ 


4.2. Spectral theory. 


Proof of Theorem^ The existence of a fundamental solution Ht is the content of Theorem 
[T^ so it suffices to show that Ht has the properties stated in Theorem Since Ht and 
dtHt are formally self-adjoint (i.e. symmetric), to show that they are self-adjoint it suffices 
to show that they are compact operators. But indeed they are, as follows from |33[ Thm 
VI.23-24] together with 

Ht,dtHt G L2(End;M X M), 

where, given a smooth section A of End, then A G L^(End; M x M) if 


^4(p, g)|||nddVolM(p)dVolM(g) < 00 . 


For t > 0, Ht is given by an H integral kernel, so is a compact operator; indeed, by (3.57), 
the index set J- of Ht G Aphg{M x M) restricted to t > 0 constant is J^(lf) = f(lf) and 
J^(rf) = f (rf), for £ the index family of H. From (3.58), these satisfy the lower bound 


inf J'(lf),infJ^(rf) > -^ -f 1 
(meaning Ht is a bounded endomorphism) and 

dVolAf(p) dVolM)?) — x^^x^k tlx dxdydydzdz, 


(4.13) 


(4.14) 


so the kernel of Ht is square integrable. Since the restriction of dtHt to a fixed time t has 
the same index set on M x M as Ht, it is also compact. 

It remains to establish (|^, i.e. that Ht{a) G T’min for every a G In fact, Ht{a) is a 
polyhomogeneous distribution with index set £(lf). This is straightforward: writing the 
expansion of iLj at x = 0 up to some order N we have 

Ht= ^ a;® log^(x)as,p(y, z, w) -f En 

(8,p)(ie(\i) 

\s\<N 
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where w = (x, y, z), and the coefficients as,p are polyhomogeneous endomorphisms on the 
manifold with boundary dM x M and is a polyhomogeneous endomorphism on M x M 
with Em = o{x^). Thus 


Ht{a) 



( ^ a;® log^(ai)as,p(z, y, w)a{w) + ENa{w))^ dVolg(w). 

(s,p)G£(lf) 

SRs<A 


(4.15) 


For example by [251 Proposition 3.20], since the x ^Em are given by a polyhomogeneous 
integral kernel, they define bounded maps of and the conormality estimates (see ( jot - 
(A.4)) follow by differentiating x~^Em- The integrals coming from the partial expansion 
terms are finite and give the expansion coefficients of Ht{a). This shows that E[t{a) € .Aphg, 
and moreover that the leading order term has no logarithmic factor. Thus, In view of Lemma 
16 it suffices to prove Ht{a) S T’max in order to conclude that E[t{a) € T’min- But indeed. 


inf (If) satisfies the lower bound (4.13), hence it follows that the lowest order term in the 
polyhomogeneous expansion (4.15) is of order at least — ^ + 1 which by (4.2) is sufficient to 
conclude Ht{a) G Because Ht is a fundamental solution of the heat equation, it follows 
that A®iLt(a) = —dtHt{a) which by the same argument is contained in since dtHt has 
the same index set as Elt for t > 0. □ 


It now follows that the fundamental solution Elt from Theorem is in fact the heat 
kernel in the following sense. 

Proposition 17. The heat kernel exp(—fA®) defined by applying the spectral theorem to 
the self-adjoint operator (A®, I?) has Schwartz kernel equal to the fundamental solution Elt 
in Theorem \l^ meaning 

a){w) = I E{t{w,w,t)a{w) dVolg{w). 

JM 

Using this we may finish the proof of Theorem 


Proof of Theorem^^ As discussed below the statement of Theorem [^ Theorem [^ itself es¬ 
tablishes essential self-adjointness of A®. It remains to prove that the spectrum is discrete, 
but this follows immediately from the spectral theorem and the fact that Ht is a compact 
operator (hence has discrete spectrum.) 

Moreover, the Weyl asymptotic formula in (1.6) follows from the standard heat kernel 
argument in 


§8.3] together with the heat trace asymptotics in Corollary 18 


□ 


Corollary 18 (Heat trace asymptotics). For each t > 0, the fundamental solution Ht in 
Theorem 12 is trace class and satisfies that F(t) := Tr Ht is a polyhomogenous conormal 
distribution on ffi"*" satisfying 


F{t) = Vol(M, g) + (^ ajf-”/2+^) + 

i=i 


(4.16) 


The proof of Corollary which uses Theorem and Melrose’s pushforward theorem, 
is deferred to Section lA.21 below. 
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4.3. Harmonic forms and Hodge theory. We begin our discussion of Hodge theory by 
pointing out that elements 7 G satisfying = 0 , admit asymptotic expansions at the 
boundary of M. Indeed, for such forms 7 , by the spectral theorem and the fact that Ht is 
the heat kernel (Corollary [T7|, we see that 


'Y = Htj = / Ht{w,w,t)'y{w)dVo\g{w). 

Jm 

By the proof of Theorem]^ specifically (4.15), we have the following. 


(4.17) 


Lemma 19. Assume that 7 G ker(A®: — > L^). Then 7 is poly homogeneous conormal 

and 7 = 0(1), i.e. is bounded in norm. 


Lemma 19 allows us to conclude that the kernel of A® is equal to the Hodge cohomology 

in ([^ . 

Lemma 20. Notation as above, 'HL^{M,g) = ker(A®: L^). 


Proof. If 7 G 77^2 {M, g) then 7 is in the maximal domains of both d and S, and so for 
smooth compactly supported /3, 

(A® 7 , / 3 )i 2 := ( 7 , A^P)l 2 = ( 7 , d 5 / 3 )L 2 + ( 7 , M^)l 2 =0 + 0 = 0 , 


so 7 G ker(A®: 

we can integrate by 

parts to obtain 


On the other hand, if 7 G ker(A®: —>■ L^), then by Lemma 


19 


0 = (AS7,7)i2 = ||d7||i2 + WhWh, 


so 7 G ni 2 {M,g). 


□ 


We can now follow the arguments in [T7|[T9] to prove Theorem]^ above. Before we begin 
we recall some facts about intersection cohomology, a cohomology theory that applies to 
stratified spaces. We do not attempt to make a full explanation of it here, but mention only 
that there is in fact a family of intersection cohomology groups for our stratified space X 
defined in ( |1.8| ) (obtained by collapsing the boundary of DM over the base +) depending 
on a function p: N —N called the ‘perversity,’ which is non-decreasing and whose values 
matter only on the codimensions of the strata of X. Here we have only one singular stratum, 
Y C X, the image of the boundary dM via the projection onto X, and its codimension is 
/ + 1, where dimZ = /. The ‘upper middle degree’ perversity m is a special example of a 
perversity, which satisfies 


m(/ + l) = 


(/ - l )/2 if / is odd, 
//2 — 1 if / is even. 


(4.18) 


The ‘lower middle perversity’m differs from rh only when / is even, in which case ni(/ +1) = 
f/2. As we will rely on the spectral sequence arguments from |17L 119) during the proof, we 
will only need to study the intersection cohomology locally, specifically on a basis of open 
sets of X. Concretely, from m, for canonical neighborhoods U = V x Ci{Z) a,s in ( |2.4[ ) 
with contractible V, we have 


From the Witt condition ( |1.4| ), we see that 
IH^iU) = IH^U) = 


if p < / - 1 - p(/ + 1 ), 

ifp> / -l-p(/ + l)- 


r m[Z) Up < f/2, 
1 { 0 } ifp>// 2 . 


(4.19) 
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regardless of the parity of /. 

Proof of Theorem^ Although Theorem describes a relationship between the Hodge co¬ 
homology and the intersection cohomology, to prove it we go through the standard route 
and use the intermediary of L^-cohomology. Thus consider the chain complex 

... ^ LlnP-\M,g) LinP{M,g) LlnP+\M,g) ^ . (4.20) 

where L^flP{M,g) is the maximal domain of the exterior derivative d, specifically 

Ll9F{M,g) = {ae L^9F{M,g) : da G L‘^nP+\M, g)}. 

Then the L^-cohomology is the quotient 

^ {dg:rjGLlnP-^M,g)} ' 

As explained in m page 6], it suffices to show that 

L2i7P(M,g)~/iJP(A;K), (4.21) 

for then the L^-cohomology is finite dimensional, which implies that the range of d (and 
thus its adjoint S) is closed. From [ITl §2.1] it then follows using the Kodaira decomposition 
theorem that 'H^ 2 iM,g) is isomorphic to L'^HP{M,g) and thus by (4.211 Theoremholds. 

Thus it suffices to prove (4.21), and for this we also follow the arguments in m pp- 5-6], 
where it is explained that it suffices to show that for canonical neighborhoods U = Vx Ci{Z) 


as in (2.4) with contractible V, the local chain complex 

LinP{U,g) 




{U,g) LinP{U,g) LlnP+\U,g) ^ ..., (4.22) 

satisfies 

L^HP{U,g)^IH^{U) (4.23) 

Here LfHP{U,g) is defined as above with U replacing M. The intersection cohomology 
groups for m are computed in (4.19), and thus we need only to analyze the groups on 
the left. To see (4.23) we use the Kiinneth formula of Zucker, [38l Corollary 2.34], whose 


assumptions are satisfied here by the fact that the exterior derivative on Z is closed on its 
maximal domain. Thus, in the notation of m Page 5], we have 


L^HP{U,g) = 0 WW{{0, l),dx^, k{p-i- f/2))) ® HP-\Z; 1 


(4.24) 


i=0 


where HTf*((0,1), dx"^, a) is the cohomology of the complex 

0^ {x‘^L^n°{{0,l),dx^)) ^ x‘^L^n\{0,l),dx^) ^0, (4.25) 

where the space on the left is the maximal domain of d on x“L^H°((0, 1), dx^). Again from 
[T^ (via [H]), WH^{{0,l),dx'^,a) = 0 if a ^ 1/2 and 1TH°((0,1), a) = M if a < 1/2 
and {0} if a > 1/2. When i = 1 , k{p — i — f/2) ^1/2 since fc > 1, so the i = 1 terms do 
not contribute. When z = 0, we have k{p — i — f/2) = k{p — f/2) which satisfies 

//2) < 1/2 if p <//2 and A:(p -//2) > 1/2 if p >//2. 

Using the Witt condition then gives 

HP{Z) Up < f/2, 

{0} ifp>//2, 

□ 


L^HP{U,g) = 


matching (4.19) and completing the proof. 
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We now discuss the proof of Theorem As the spaces in the theorem are incomplete 
cusp edge spaces in a neighborhood of the divisor by [2S], our results would apply to these 
spaces, if not for the fact that moduli spaces such as these have interior orbifold points. 
This is not a problem, since, as in |20] we may lift to a finite cover with no such points. 
One can then work on the space of functions which near each orbifold point 

are smooth when lifted to a local finite cover resolving the singularity. Constructing a heat 
kernel on the lift and averaging over the group action then gives a fundamental solution to 
the heat kernel downstairs which has all the desired properties. We leave the details of this 
simple extension to the reader. 


Appendix A. Manifolds with corners 


In this section we recall some of the facts about distributions on manifolds with corners 
(mwc’s) used in this paper. This material is due largely to Melrose, and the reader is referred 
to his book m for more details. See also [IB] . 

The objects considered here, for example the ice-metrics (2.9), have polyhomogeneous 
regularity, which we define now. The sheaf of polyhomogeneous conormal (or polyho¬ 
mogeneous, or simply phg) functions Aphg(A') is defined as follows. First, an index set £ on 
a manifold with corners X is an association to each boundary hypersurface H oi X a set 
£{ H ) C C X N satisfying that the subset 
{{z,p) G £{H) : Re 2 < c} is finite for all c G M. ^ ^ 


Given an index set £, for a boundary face F = for boundary hypersurfaces Fli, 

define the subset £{F) C x W by (z,p) = {zi,..., zs,pi, ■ ■ ■ ,ps) G £{F) if and only if 
{zi,pi) G £{Hi). We define the Frechet space Api^g{X) as follows. We write u G ^pjjg(X) if 
and only if for each boundary face F = writing pi for a boundary defining function 

oi Hi, u satishes 


u ^ o.z,pP^ log^ p where 

(z,p)&S{F) 

=w pt > log^ ^=n Pi ’ 

i=l 

and the symbol ~ means that 

Ej^ = u- y] log^ p, 

Re zi<N Mi 


(A.2) 


(A.3) 


where Em is a smooth function on the interior of X which is 0{\p\^), where \p\ = {p\ -I- • • • -I- 
pIY^'^. Moreover, Em is conormal, meaning that if Vb = Vb{X) denotes the set of smooth 
vector helds on X that are tangent to all boundary hypersurfaces, then 

\p\-^V^E c L°°. (A.4) 


Note that if a phg function u vanishes to infinite order at H, then u is polyhomogeneous 
with index set £ satisfying £{H) = 0. 


Lemma 21. Let X denote a mwc, Xi{X) = boundary hypersurfaees, and for 

each z G I, let pi denote a boundary defining function of Hi. Given a smooth vector bundle 
E —> X, if Ki are polyhomogeneous sections on Hi, then provided 

PT^jlHinHj = Pj' 


(A.5) 
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there exists a polyhomogeneous conormal distribution K on X satsifying 

p1'K\Hi = Hi (A.6) 

Assume moreover that at a particular boundary hypersurface which we take to be Hi, that 
we are given an index set Fi C C x N and polyhomogeneous sections bj^p G Hi). 

Then given functions Ki on Hi, i 1, there exists a distribution K satisfying ( |A.6[ ) fori 1 
and such that 

^ P? logP(pi)&s.P (A.7) 

S ,PG Fi 

provided (A.5) holds for i,j^l and furthermore for i 1 

^i-PT Pi (^-8) 

s,peFi 

Remark 22. (1) Note the converse; if K = a for some positive function a near 

Hi n Hj then setting = ki for I = i,j, we have p’^f Hi = Pi'Kj on Hi n Hj. 

(2) The matching condition (A.5) implies further matching conditions on multifold in¬ 
tersections, e.g. it implies that 

Pi'p"/ Ki = pl'p^Kj = pI'^p''/ Ki on Hi n Hj n Hi. 


(3) The second matching condition (A.8) merely says that the desired data on a bhs Hi 
has the same asymptotic expansion at Hi as the the desired distribution restricted to Hi. 


Proof. Denote the number of boundary hypersurfaces of X by m = |Af |. There is a number 
<5 and boundary defining functions pi such that the set {pi < 5} is diffeomorphic as mwc’s to 
Hi X [0,(5). Without loss of generality we take J = 1. Following the remark, for a collection 
of bhs’ Hii,..., Hi ^, the distrubution 

Kii...ip = ( P |pij=--.=pip=0 

iy^ik 

is well-defined independently of the choice of 4- G to}. 

Let x{^) be a cutoff function with x = 1 for a; < 1/3 and x = 0 for a; > 2e/3. For the 
distribution K we may take 

Ttl 

K = n x(Pa)p7''')- 

P=1 l<ii<'"<ip<m iG{ii,...,ip} 

For example if to = 2 then 

K = x{Pi)pT^Ki + x{P2)P2‘"'k2 - x{Pi)x{P2)pT^pT'' 

Note that each term in the sum defining K defines a polyhomogeneous conormal distribution 
on all of X, as the distribution Kii,,,i is defined on a neighborhood of Hi^ D ■ ■ ■ O Hi off 
which the product i } l^(Pi) vanishes. 

Letting be the term corresponding term in the definition of K, note that if z ^ 

{ii, ...,ip} then p‘l'Aii,,,i^ = pf |p,=o- Fixing i, multiplying by p^'K and restricting 

to Pi = 0 gives 

m— 1 

prxu=o = Az. + ^(-i)*’-vr( E E ^u...zp+Ju=o 

p l<ii<-"<2p<m l<ii<---<2p+i<m 

i^{ii...Zp} iG{ii...ip+i} 


= Ki, 
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which establishes (A.6|. 


We now prove the final statement of the lemma. Let x be the cutoff function defined 
above. First, we claim that under the stated assumptions there exists a distribution K' 


supported in {pi < 1} satisfying both (A.7) (with K replaced by K') and that 


pTK'\h, = x{Pi)Ki (A.9) 

for each i ^1. To see this, take any distribution K" supported in {pi < 1} satisfying ( A.7| ), 
and note that at := p^'K"\Hi —x{Pi)i^i = By the support condition, the distribution 

K' = K" — x {Pi)^i is defined globally, has the same asymptotic expansion at Hi as AT", 
and satisfies (A.9). This K' will play the role of xipi)p ~^^from the previous paragraph. 
Concretely, for 1 < < • • • < ip < m, let Then we 

may take 


K 


E(-i) 


p-i 


x{pj)p] 










+K'+j2{-ir-^ 

p^l 


l<ii<---<ip<m je{ii,...,ip} 


( n X(Pj)p/0- 


Again, for example if m = 2 then 

K = K' + x{P2)pT^i^2 - {p'2^K')\h^p~‘''^x{P2)- 
The given expression for K can be directly checked to satisfy (A.6)- 


□ 


A.l. Melrose’s pushforward theorem. Given a map /3: X —Y between manifolds 
with corners, if At (•) with • = X,Y denotes the space of boundary hypersurfaces, then /3 
is a b-map if it is smooth and if for each H € Xi{Y) with pn a boundary defining function 
for H then 

■Ph' 


P PH = aRH'(^M(x)PH' Pm 


where a S C°°{X) is non-vanishing and N is the number of boundary hypersurfaces of Y 
and the are non-negative integers. This means foremost that pH pulls back to a 

smooth function, and the numbers e{H',H) simply keep track of the order of vanishing of 
P* Ph at each face of X. The function 


e-. M{X) y. M{Y) —^ No (A.IO) 

is the exponent matrix of /3, and e{H', H) > 0 means H' maps into H via /3. 

If a b-map has a few additional properties then it pushes forward polyhomogeneous dis¬ 
tributions (more accurately, densities) to polyhomogeneous distributions and their index 
sets change in a way dictated by the exponent matrix. Note that it follows from the defi¬ 
nition of a b-map that every boundary face A of A (meaning an intersection of boundary 
hypersurfaces), can be associated to a face P{F) of Y defined to be the unique face with 
P{x) G P{F)° for every x G F° . A b-map /3: X —F is a b-fibration if: 

• P does not increase the codimension of faces, i.e. for each boundary face F of X, 
the associated face P{F) in Y satisfies that codim(A) < codim(/3(N)). 

• Restricted to the interior of any face P: F° — > {P{F))° is a fibration of open 
manifolds in the standard sense. 

According to a theorem of Melrose m which we state below, a b-fibration pushes forward 
phg densities to phg densities in a manner we describe now. First, on a manifold with corners 
we choose a non-vanishing b-density p, meaning a section of |A|"(^r*A), the density bundle 
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of the b-cotangent bundle. The b-tangent bundle ^TX is the bundle whose smooth sections 
are Vb, the vector fields tangent to the boundary. The bundle ’>T*X is the dual bundle of 
^TX, and near a face F = where pi are bdf’s and y are coordinates on F then, the 

sections of ^T*X take the form 

It follows that near any intersection F = Dj^jFlj^ of boundary hypersurfaces for J C I 
where I indexes Xi{X) (i.e. any face of X) that a non-vanishing b-density takes the form 


M = 


dyUn(^jdp. 


rijeJ Pi 


(A.ll) 


for some smooth non-vanishing function a on X. A polyhomogeneous b-density u S 
Apijg(Ar) (g) \A\‘^{^T*X) can be written as fp for a phg function / and the index set of 
u is by definition the index set of /. 


Theorem 23 (Melrose ISHl)- Letu € Mpj^g(A')(g)|A|”(^T*Ar) be a polyhomogeneous b-density 
on X with index set £, let f: X —> Y be a b-fibration with exponent matrix e, and define 
the pushforward f^,u to be the distribution on smooth funetions v G o.cting by 

{f*u,v)Y = {u,f*v)x- Then provided for each F[ G A4(X) we have 

e{H, H') = 0 Vi/' e M{Y) £{H) > 0, (A.12) 

then f^u G Vlphg(T) g) |A|”(*’r*y) where 

with the (extended) union taken over FI' with e(ii',ii) > 0. 


The extended union, defined in EH, contains the standard union and possibly more log 
terms. 


A. 2. Heat trace asymptotics. We now use Theorem to prove the heat trace formula 
in Corollary [T8| above. The heat trace is equal to 

Tr(e-‘'^)=/ Ht{w,w)dVo\g = a4{L*Ht)dVo\), (A.13) 

Jm 

where l: M x [0, oo) — y M x M x [0,oo) is the diagonal inclusion and a: M x [0, oo) —>■ 
[0,oo) is the projection onto the right factor. The natural space here on which to consider 
Fit is and thus to evalute this pushforward we must see how <j and t act on the 

natural blown up spaces. The following may be easily verified. 

(1) The closure (Ai^ggj)A := cl(6(M° x (0, oo))) is a manifold with corners with 4 bound¬ 
ary hypersurfaces, sf, ffi'^, ff'^, tf'^, equal to the intersection of cl(6(M° x (0, oo))) with 
rf n If, ffi, ff, and tf, respectively 

(2) The map a extends form the interior M° x (0,oo) to a b-fibration cr: — >■ 

[0, oo) with exponent matrix 

e,(sf) = 0, e<,(fr/) = 2, e,(ff'') = 2k, e<,(tf‘^) = 2. 


To apply the pushforward theorem, we note that the volume density 
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is equal on to 


fi = a {psfPs,<‘Pgd)’‘f+^^lo, 


where po is a non-vanishing b-density on Thus {l*H)p is phg on with 

index family £‘^ satisfying 

inf£’‘^(sf) = 3, inf £’'^(ff/) = -6, 

inf £‘^(fr‘^) = k{f -n) + l, = {-n, -n -f 1,... }. 

Note that Tre“*^Y = Ht)p)- The integrability condition (A.12) must be checked 

only for sf and thus holds by Theorem |12[ and we apply the pushforward theorem to obtain 
that Tr is polyhomogeneous with index set 

{(Ci/ 2 ,pi): (Cl,Pi) e £:‘^(fri")} u : (C2,P2) g 

DKCsAps): (C3,P3)G ^(tf")}- 

In particular, 

/ 

Fit) = +cot-(^+i)/2+i/(2fc) +0(i-('>+i)/2+i/(2fe)+.)^ 

3=0 

for some e > 0. As discussed in EZl Section 3.3], the heat kernel in fact lies in an even 
calculus and thus the terms for odd j in this sum are equal to 0, giving the trace formula 
(4.16). The fact that the leading order term is the volume is standard. 

Appendix B. Triple Space 

We will now analyze composition properties for “Volterra” type convolution operators 
as described in (3.70). To do so, following [SI IS], we construct a “triple space,” which 
we denote by MAj]which is designed specifically to accomodate the process of composing 
operators which have the structure of the error terms in (3.71). The structure of our triple 
space is analogous to that constructed by Grieser and Hunsicker in with slightly different 
homogeneities and with the added complication that there are time variables involved. 
Note that, given Ai, i = 1, 2, we want is to make sense of the integral 

f f Ai{w,w\t')A 2 {w',w,t — t') dVolg{w')dt'. 

JmJo 


Define the wedge 

W := {t-t' >0} C 

and define the left, center, and right projections 
ttl : M X M X M X W 
(w, w', w, t, t') 
nc ■ M X M X M X W 
(w, w', w, t, t') 

TTfi : M X M X M X W 
(w, w', w, t, t') 


p+ 
h' j 


(B.l) 

(B.2) 


M X M X [0, oo)t 
M X M X [0, oo)t 

iw,w,t) 

M X M X [0, oo)i 

(w', w,t ~ t')- 


(B.3) 


Then, formally, the integral in (B.l) says that the integral kernel of A 1 A 2 (as an operator 
acting by convolution in time) is 

{AiA2)iw,w,t) = (7rc)*(A-4i)(AA), (B.4) 
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where (ttc), denotes the pushforward, i.e. the integral along the fibers of ttc (which, by the 
way we have set up the problem, requires the choice of a metric on the fibers which we come 
to shortly.) Analysis of (B.4| becomes tractable if the space x is blown up so that 
the pushforward theorem described in Section [A. 1| applies. 

Note that x bb is a manifold with corners with 5 boundary hypersurfaces 

L = {a; = 0}, C = {a:' = 0}, R = {S = 0} 

tb’i = {t' = 0}, tb ’2 = {t — t' = 0} 

It is easy to check that, in the language of Appendix [A] the maps tt, with • S {L, C, R} are 
b-maps from x bk to x [0, oo)t and the exponent matrices are also easy to compute. 




• = L,.' = If 

• = C, •' = rf 
• = tb^, •' = tb 

otherwise 




(•>•') = 




• = C, •' = If 

• = R, •' = rf 
• = tb2, •' = tb 

otherwise 


• = L,.' = If 

• = R, •' = rf 
otherwise 


(B.5) 


(B.6) 


We blow up X bk to form a space (3: M? 


heat 


X bk in a sequence of steps as follows. 
First, consider the three pullbacks of the submanifold Bq = {x = x,y = y,t = Q} G 

(B.7) 


X [0, 00 )^ defined in (3.1) 

TT^h^o), ^r\B^). 

These three sets intersect pair-wise in the triple intersection: 

n TTch^o) = Trc\Bo) n 7r^i(So) = TT^h^o) n ^-r\Bo) = 5, 

where 

S = {x = x' = x = t = t'=y — y' = y' — y = Q}. 

We blow up the set S, with appropriate homogeneities, specifically letting 

-^heat.O ~ ^ bF;5]i„/iomj 


(B.8) 

(B.9) 

(B.IO) 


with t x^ 


ix') 


/\2 


\y - y'\ ^ \y' - y\ , and let 


I • -^^heat.o 


X bF denote 


the blowdown map. Call the introduced boundary hypersurface ffi . Near to ffi , we have 
polar coordinates 

pn = (t + + (x')^ ^-x"^ + \y-y'y + \y' - y^) , 

t' t — t' X x' X y — y' y' — y 
Pn ’ Pn Pn 


2 ’ 
Pn 


Pn 


(B.ll) 


=: 


Yx' J Yx 


Pn Pn 
,y'y-yn(l>y>-y), along with y',z,z',z. 


v-y ’ -y/ 

The asymmetry of the y,y',y in the coordinates is spurious in the sense that if one defines 
^’^-y = {y - y)/pn, then any two of the can be used in by redefining pn 


using e.g. \y — y'^ and \y — y^ (and then using Either set of coordinates is 

defined in a colar neighborhood of ffi'^. 

We then blow up the closures of the lifts 

5. :=cl((T.o^o)-'(So)\ffi^), 
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i.e. the rest of the lifts of the Bq via the three projections, where • S {L, C, i?}. These are 
disjoint subsets and we blow them up in any order, setting 


-^heat,l “ [-^heat.oi ^•=L,C,RS»\inhom 

with the appropriate homogeneties, e.g. for Sl have t' ~ 
we have a blowdown map 

Pi: 1 —> M X M xM xW. 


(B.12) 

(x')^ ~ Ij/- 2/T- Again, 


(B.13) 


The new faces we call ff* with • S {T, C, i?}. Coordinates at ff^" can be determined as 
follows. Note that Sl is given in the coordinates (B.lll by = pQ = pQ, = (l>y_y' = 0, 
and that in a neighborhood of away from ffi'^, ~ t'. Thus, to match homogeneities 
with the blowups of the double space, we want to blow this up so that the following give 
polar coordinates near the intersection of and ffi'^: 


y-y 


(^pL)2^ pL^ pL^ pL 

-■ along with 


(B.14) 


with functions as in (B.ll). It is also possible to use simpler projective coordinates, as we 
will see below. Coordinates near Sf can be derived similarly by switching t/fp with pQ-t' 
and with (/)-. The situation at fff is slightly different since, writing = (j)^, the 

pullback of (j)'^ on ^ via ttc vanishes at (/)p = 0 = and thus Sc is codimension 1 

higher than 5, for • = L, i?. 

Here we blow up so that the following give coordinates 

= { 4 '? + + \ 4 y-y\^) ^ ) 

( 4v 47-t' 4Q 4~ 4y-y 


irf = 


V(p^)2’ pC 


(B.15) 


= : {(t)t,,4i-t'^4x^4'iAy-y), ^^OTLgwiihy' ,z,z' ,Z,Pn,4l2'^4y'-y- 


Lemma 24. With terminology as in Appendix \A.l\ the maps tt, extend from the interior 
to b-maps 

■ Afheat.l ^ Afjjeat,! 

for • G {L, C, i?} with exponent matrices e??, satisfying 

es.(fr?,fri) = l,es.(fr*',fri) = 5.,.., (iff, if) = l,e^^(fff,rf) = 1, 


(B.16) 


= l,esi(fff,rf) = 1, eSfl(ffi', tb) > 0, (Iff, tb) > 0, 


(B.17) 


where 5,,,i =1 if • = •' and zero otherwise. When • € {L, C, R, tb^, tb^}, i.e. when it is 
the pullback of a boundary hypersurface of MxMxMxW via the blowdown map, then 
the exponent matrix satisfies (B.5) with tt replacing tt. 

Moreover, ttc is a b-fibration in the sense of Appendix\A . l\ 


Remark 25. The significance of the inequalites in (B.17) involving tb is that all the dis¬ 
tributions under consideration vanish to infinite order at tb, and thus the pullbacks of these 
distributions via will vanish to infinite order at S'f, and the same for tt^. 
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Proof. We verify the lemma for for nc and leave the other nearly identical calculations to 
the reader. That ttc extends to a b-map follows easily by writing the pulling back the 
coordinates in (3.3) and writing them in terms of those in (B.ll). In particular, note that 
the pullback 

^cP = PnP^, (B.18) 


so the exponent matrix claim holds. The rest of the definitions of b-hbration are easy to 
check. □ 


^ c 

Remark 26. The extended map ttl is not a b-fibration as it maps the interior ofSf to the 

interior of the face tb H If due to the fact that t = 0 on W implies that t' — 0 also, thus the 

map increases the codimension of a face. The same holds for tt/j, i.e. 7f/j(ff^) C tbn rf. 


Next we must blow up the lifts of Bi in (3.4). Since by (B.17), tt, only maps If* to ffi 
if • = any of the pair-wise intersections is again equal to the triple intersection 

5 ' = n nc\B^) = n n 

Indeed, each is a subset of ff^, and in the polar coordinates defined on the interior of ff^, 


using the definition of Bi in (3.4) 


5' = {p = = 0, = cjff, = cfq,4ry_y, = cjp^_~ = 0} , (B.I9) 

with no restrictions on y',z,z',z. We form a space S']inhom with appropriate 

homogeneities. To understand this space, note hrst that near S' we can use projective 
coordinates on ffi , concretely we can take for example a; to be a boundary defining function 
of ffi n and coordinates x,t' jx?', {t — t')fx'^,xlx,x'fx, {y — y')/x, {y' — ^/x to replace the 


polar coordinates in (B.ll). Then the homogeneities are determined by those in the ff 


blowdown of the double space, and one has coordinates 




X — X, 


Pn 


^ 1 ^x—xi — XI — y) 




l/2(k-l) 


t' 


t-t' 


X — X 


X -X y-y y -y 


n ’ ~n-2(k-l) ’ ~ 9 _ 2 (fc-l) ~_(fc-l) ’ ~_(fe-l) ~_(/c-l) ’ ~_(fe-l) 

Pn x^Pft Pn ^Pn ^Pn ^Pn ^Pn 


along with y, z, z', z. 

(B.20) 


One can also take coordinates in which x, x', x are permuted, and the same with y, y', y. 

We let ff'^ denote the introduced boundary hypersurface. 

The lifts of the n~^{Bi) minus their intersections now have disjoint closures. For example, 
we have 

n ffl^ \ff^ = {f,, = = 0} 

while 

n ffi^ \ ff^ = R = = 0}, 

whereat = and and for ttr we have = V't-t' = 0, tpy' = 0; 

since \4>\ = 1, these sets are disjoint. Furthermore, the pullbacks satisfy that 

n-\Bi)nSl =6, y, 
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for € {R,C,L}, and each intersection is straightforward to write down, e.g. with coor¬ 
dinates as in (B.15), 


s^' = {p^' = - 0? = 0 ,= 0}. 


y-y 


We will blow up first the n ff^ and then the n~^{Bi) n ff* with for • G {L, C, R}. 

In the interior of ff* with • G {L, R} the blow ups of the pullbacks of Bi are particularly 
easy to understand as there we can just pullback the projective coordinates in (3.91 and use 


these together with the other uneffected coordinates to obtain projective coordinates e.g. 
near ttl o /3q^(Bi) H fff valid near the interior of the introduced boundary hypersurface. 


s — 1 


(x') 




(x') 


f\k 


/ y 

V = 


x'Y 


T' = 


t' 


{x') 


,/\2k ' 


(B.21) 


together with w, t on the introduced boundary hypersurface. In the interior of ff^, one 
needs only to remember that the vanishing of the pullback of the ij)t coordinate implies the 
vanishing of both Yt' and Yt-t' ■ One can use a? as a boundary defining function and then two 
projective time coordinates T' = i! jx^^ and T = {t — t')/x^^- In the interior of ff^ but away 
from ff*^, we want the same homogeneities, but now the pullback of x' in the interior of ff*^ 
is proportional to and in the interior of the introduced blow up we will have coordinates 
as in (B.21) with all the functions replaced by their </) couterparts, e.g. x' replaced by Yx' 


and ^ 7 ^ replaced by ipy-y'/Yx 


We focus at the intersection ffi'^ n ffi*, first with • = C. Near Sc, we can simplify things 


slightly by using projective coordinates, derived from (B.15) by noting that (f>~ is non-zero 


at ffi'' n ffi* n Sc and can thus be used as a boundary defining function. Specifically, take 


x = x = T=^, y = ^, 

X X^ X 


together with the other (non-polar) coordinates in (B.15). Blowing up to introduce a face 
ff'^’*^, have 

V = {T + pT~'^'^ + {X- 1)^ + |;y|2)i/(2(fe-i))^ 

but it follows that iSi n ffi*^ intersects ff'^’^ at vk = (0,1, 0, 0) and thus p'~' can be used as a 
boundary defining function. Again working near Sc we can take pn as a boundary defining 
function for ff'^’'^ and use projective coordinates pn,T/pn*^~^\ (df—l)/pn“^,(V/pn“^. Using 


these we blow up H ffi*^ with 


V = { 


T 


2 (/c-l) 

Pn 




-11 , (-r -1)^ 


k-1 

Pn 


l^_!l)l/(2(fc-l))^ 


k 

Pn 


^ = (7^ 


T 


^ x-i _ y 

(Ppn)2(/C-1)’ r ’ 




X — X 


y-y 


x‘^{VpnY^^~'^^ ’ Vpn ’ x{VpnY-'^ ’ x{Vpn) 


k-l ' 


and this is the final blow up of Sc- The blowups for Sl^Sr are similar and left to the 
reader. 
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Proposition 27 (Incomplete cusp edge heat triple space). The above construction yields a 
space and blowdown map 

(B.22) 

heat —^ ^heat ,2 exponent 


j3: —y M xMx M xW, 


such that the maps tt, from (B.16) extend to b-maps tt, : Mf 


(B.23) 


matrix satisfying (B.5|, (B.17) (with tt and tt replaced by tt), and 

= 1, = <5.,.- 

= = 1 

e^,(ff^,rf) = = 1 

Moreover, apart from components o/e^p, •') and •') with* € {L,Ii}, all other 

components are zero. 

Again, nc is a b-fibration. 

Proof Again, we focus on ttc- To check that nc extends to a 6 — map, we pull back the 
polar coordinates p, (j),y, z,z, from (3.7) defined at ff in First, we compute 

+ (s -1)^ + i\y- y\/x)^) ) 


= (t + + {X- If + 13^1") 


s l/2(fc-l) 


= Vp, 


C) 


and then note that 


x X — X 


y-y 


t X X — X 


y-y 


□ 


,a;2(Ppn)2('=-i) ’ pp^ ’ xfP ’ xifP 

= A). 

This establishes both claims for ttc- The R,L case are left to the reader. 

Proposition 28. For i = 1,2, let Ai G “4phg(-^heat 2 ) index sets £i satisfying 

£i{^i) > —3 — b— kf, £i(ff) > —kn — 2k, 5i(lf) = 5i(tb) = 0 and £i{ri) satisfying (3.57). 
Then A 3 := ff* Ai(w,w',t')A 2 (w',w,t — t') dVol^,,' dt' lies in f where for any 

e > 0 , 

inf £ 3 (ffi) > inf £’i(ffi) + inf faCffi) + 3 + 6 + fc/ - e, 
inf £ 3 (ff) > inf £(i(ff) + inf £( 2 (ff) + fcn + 2 A: — e. ^ ^ 


Remark 29. The constants kn + 2k and 3 + b + kf in (B.24) should be interpreted, for in¬ 
stance in the case o/ff, as saying that the (Volterra type) composition of two operators given 
by Schwartz kernels as in the theorem has Schwartz kernel whose leading order asymptotic 
behavior at ff increases relative to the rate —kn — 2k, in particular if both the composed 
operators grow like —kn — 2k then so does the composition. These are, incidentally, the 
exact rates of blowup of the heat kernel times t~^ at the faces ff and ffi, and furthermore the 
fact that the errors t~^Q in (3.71) vanish one order faster than t~^H means, as described 
above, that taking powers makes them vanish at increasing rates at both ff and ffi. 
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Proof. We write A 3 as the pushforward of a b-density and then apply the Pushforward 
Theorem from Section A.l First we define a non-vanishing b-density on Mx M xMxW 
as follows. We let 1 / be a non-vanishing b-density on M satisfying ix = a \| for a smooth 
nonvanishing function a near the boundary, and consider 


dt'dt 


Hq = V V V 


— t') 

where ix',lx are equal to ix in the primed and tilded coordinates, respectively. Since the 
blowdown map (3 from (B.22) is a b-map, is a b-density on and one checks that 

pyo = Gy, (B.25) 

for a non-vanishing b-density y on and G G G°°{M^^y satisfying that for some 

non-vanishing smooth function G', 

G = G'pg^cp^^n{pffL pgc^ 

X {psn.Lpgn,cpsn.uy+^'>^+>^-yfn^+^^^~^\ 

Then we can write the desired pushforward as a pushforward of a b-density, specifically 


^3 \j\) = (ttc)* (TT^AiTrJjAa • ((t'A)(i “ t'))F{w')po) 


(ttc)* (wlAiTr*j^A 2 ■ P*{{t'/t){t - t')F{w')po) 


(B.26) 


where F is the function defined by dVolg = Fix and in particular 

F = a 


where a is a non-vanishing polyhomogeneous function on M, and v are the pullbacks 
of the density v above to the left and right spacial factors oi M x M x K+. To find the 
asymptotics of A 3 itself we must compute the asymptotics of the densities on the left hand 
side of (B.26); Letting P 2 again denote the blowdown map ^ ^ M x M x [0,oo) in 


(3.8), we check that 


dt.^ 




where p 2 is a non-vanishing b-density on Thus from (B.26), if the distribution 


{'Kc)*y*LAiyA 2 ■ /t){t — t')F{w')po)) is polyhomogeneous with index set £3 then A 3 

is phg with index set £3 satisfying 


£:3(ffi) = 4(ffi) - b, £3iS) = (ff) -{kb + k- 1), 


(B.27) 


and £ 3 (») = f 3 (») otherwise. 

Thus the index family of /t){t — t')F{w')) must be determined. To 

determine iffAi, we see that, at a bhs F[ of the index set of if’fAi is simply the 

index set of Ai at the bhs H' of 2 which H is incident. Thus from our work above 
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we see that has index set satisfying 


fi(tb;) 


£i(L) =£i(lf) = 0 
£i(C) = = £i(rf) 

= £i(fri^) = = £i(tb) = 0 

£i(fri^) = = £i(fr^’^) = £i(fri) 

f = Si{S^) = £i{S) 

£i(R) =Z, 


the last line coming from the fact that is independent of x, in particular is smooth 

up to R. The index set £2 oi 7 f^A 2 has the same expression in terms of £2 but with all 
‘R’s switched with ‘L’s, all If’s with rf’s, and all I’s with 2’s (except of course for the 1 in 
the subscript of ffi). For example, (c.f. the second line above) £ 2 (C),f 2 (ffi^),'? 2 (ff^) are 
all equal to £ 2 ( 11 ), which is assumed to be 0. If we define the operation £i 0 £2 on index 
sets to denote the index set whose elements are sums of elements of the two index sets, It 
follows that 7 rJjA 2 is polyhomogeneous with index set T satisfying 

{J-(C),J-(L),^(ffi^),^(ff^),J-(tb'i),^(tb'),J-(ffi'"),-^(ff^’''),-^(ff'')} = 0 
j-(ffi^) = £:i(ffi)©£:2(ffi), J'(ff'')=£i(ff)©£:2(ff), J'(R) = £2(rf) 

/■R OQ\ 

J-(ffi«)=£:i(rf)©£: 2 (ffi), J-(ff^) = £:i(rf)©£ 2 (ff) 

j-(ff^’^) = £'i(ffi) © £: 2 (ff), ^(ff^’^) = 5i(fr) © £ 2 (ffi). 


Now we compute the asymptotics of the term /t){t — t'))F{'w')fio) = — 


t'))F{w'))G^Q with G in (B.25). First, write P*— t'))F{w')) = • 7 rJ^(i) 7 rJj(t) 7 r|j(F") where 
F is thought of as a function of the left factor of M x M x [0, 00 ). Recalling p, p from (3.3) 


and (3.71, respectively, and letting a denote a polyhomogeneous function which is smooth 
and non-vanishing up to boundary hypersurfaces • for which F{u) 7 ^ 0 (and whose value 
will change from line to line), we compute 




= a 


(Pff 1 Pff 1 Pffn ■« ) ^ (Pffn Pff 1 n A Pffi ) 


2k 


(Pff© Pffi C Pffn,HPffn, 1 . )2 (pjjn n,C PgC )2^= 

X (Pff© Pffi« Pffn A ) ^ (PfT^ Pff©.« 

X (PCPSi ^PffiPffinPffi«PffnAPffnpgn,i?pgR 
= a (Pff©)^(Pff©APffI-)^''(pff©Pff^R)^(pffnpg^n,RpgR)^'' 
X (PCPffi i- Pffi Pff© Pffi H Pffn A Pffn Pffn ,i? PgR ) 

= a (Pff©Pffi«)^(Pff'^PffnAPffn,RpgR)^'" 

X (PffiHPffin Pffn pgn,Rpgn,I. pgR) 
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Putting this all together, we see that 7r^y4i7r^y42 • /t){t — t'j)F{w')no) is polyhomo- 

geneous with index set T 

0 £2(ffi) + (3 + fc/+ 2&), 
r(ffn) = 0 0 (1 0 2fc 0 fc/ 0 2 kb), 

J=(R)=£ 2 (rf) 

= £i(rf) 0 £: 2 (fFi) 0 (3 0 fc/ 0 6), (B.29) 

= fi(rf) 0 £:2(ff) 0 (1 0 2fc 0 fc/ 0 kb) 

0 £ 2 ( 3 ) 0 (1 0 2fc 0 fc/ 0 (A: 0 l)b), 

= £i(ff) 0 £2(ffi) 0 (1 0 2fc 0 fc/ 0 (A: 0 l)b), 


and J^(*) = 0 for all other values of •. 

Now we apply Theorem 23 to analyze (0c)» ~ from 

(B.26). To check that the conditions of the theorem hold, we first recall that 0c is 
a b-fibration. Also, note that e^^{C,H') = (tb^, i?') = 0 for all 

H' e '^(-^heat 2 )i must check the inte grability condition there, but by be¬ 

low ( |B.29 ) we have F{C) = ^(tb^) = J^(tb 2 ) = 0, so the integrability condition holds. 
Thus A^{-k’^)*i' is phg on ^ ''^33 index set £3 satisfying 


£'(lf) = J(L) D D = 0 

£'(rf) = J(R) U J(ffi^) D 

£'(fr 3 ) = D D D J(ff^’-^) (B.30) 

£'(ff) = F(ff^) U D ^(3*^) = F(3^) 

£'(tb) = 0 , 


where we used from below (B.291 that various bhs’s have infinite order vanishing. From this 
we see that the bounds in Proposition |28| hold, in particular that for any e > 0, 


inf (ffi) > inf (ff 1 ) + inf ^2 (ffi) + 3 + fc/ + 26 - e, 
4(ff) = inf £’i(fr) + inf £’ 2 (ff) + 1 + A: + fcn - e, 


and thus by (B.27) the acutal index set ^3 of A 3 satisfies (B.24|, and the proof is complete. 

□ 
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